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ABSTRACT 

This  study  simulated  the  operation  of  the  T-46A  trainer 
aircraft  In  the  Undergraduate  Pilot  Training  (UPT) 
environment  In  order  to  estimate  aircraft  availability  and 
sortie  generation  rate.  The  simulation  model  Is  based  on 
current  T-37  aircraft  UPT  operations  and  uses  estimates  of 
the  reliability  and  maintainability  of  the  T-46A. 

Regression  analysis  techniques  were  used  to  estimate  the 
functional  relationship  between  the  Independent  variables  and 
the  response  variables.  After  Initial  screening,  only  two 
factors  were  Included  as  the  Independent  variables.  These 
were  mean  time  between  failures  (MTBF)  and  mean  time  to 
repair  (MTTR).  A  central  composite  design  was  used  to  gather 
the  data  needed  to  perform  the  regression. 

The  results  of  the  regression  analysis  Indicated  that 
for  both  aircraft  availability  and  sortie  generation  rate,  a 
second  order  regression  equation  in  terms  of  only  the  MTBF 
factor  provided  the  best  fit.  As  was  expected,  an  increase 
in  MTBF,  meaning  the  aircraft  is  more  reliable,  results  in  an 
Increase  In  both  aircraft  availability  and  sortie  generation 
rate.  Estimates  and  confidence  intervals  for  aircraft 
availability  and  sortie  generation  rate  were  determined. 
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A  SIMULATION  MODEL  OF  THE  T-46A  AIRCRAFT 
FOR  AVAILABILITY  AND  SORTIE  PROJECTIONS 


I.  Introduction 


Background 

One  of  Air  Training  Command's  (ATC's)  primary  functions 
Is  to  provide  flight  training  to  undergraduate  aircrew 
candidates.  Currently,  ATC  uses  the  Cessna  T-37B  aircraft  as 
the  primary  trainer  to  perform  this  function.  Unfortunately, 
"(t]he  T-37  Is  a  1950's  technology  trainer  that  Is  becoming 
Increasingly  more  costly  to  operate  due  to  high  fuel 
consumption,  price  escalation  of  fuel  and  parts,  availability 
of  parts,  and  Increases  in  manpower  costs"  (4:1).  In 
addition,  the  T-37  does  not  have  a  pressurized  cockpit  which 
limits  the  airspace  In  which  the  T-37  is  able  to  perform  its 
training  maneuvers.  Finally,  the  number  of  T-37's  available 
is  less  than  the  number  projected  to  be  needed  In  future 
years. 

Because  of  these  deficiencies,  the  Air  Force  has  made 
the  decision  to  buy  a  new  primary  trainer  for  flight 
training.  Fairchild  Republic  has  been  given  the  contract  to 


produce  this  new  primary  trainer  called  the  T-46A.  The 
objectives  of  the  T-46A  are  to  overcome  the  "...operational 


deficiencies  of  the  T-37B,  to  realize  operational  and  support 
cost  savings  through  the  use  of  modern  airframe  and  engine 
technology,  and  to  provide  ATC  an  adequate  number  of 
airframes  to  meet  flying  hour  requirements  beyond  FY  87" 
(4:3). 

Fairchild  Republic  will  deliver  the  first  two  production 
aircraft  to  the  Air  Force  In  April  1986  (7:109).  The  Air 
Force  Operational  Test  and  Evaluation  Center  (AFOTEC)  will 
use  these  two  aircraft  to  perform  Initial  operational  test 
and  evaluation.  However,  before  these  planes  are  delivered, 
the  Aircraft  Logistics  Analysis  Branch  In  AFOTEC 
( AF0TEC/LG4 A )  Is  Interested  In  obtaining  a  simulation  model 
which  will  estimate  certain  performance  characteristics  of 
the  T-46A.  Specifically,  AFOTEC  is  Interested  in  predicting 
aircraft  availability  and  the  sortie  generation  rate  for  the 
T-46A  aircraft.  Aircraft  availability  Is  the  percent  of  time 
that  an  aircraft  Is  capable  of  performing  all  of  its  assigned 
missions.  The  sortie  generation  rate  Is  the  average  number 
of  sorties  (flights)  an  aircraft  is  capable  of  performing  In 
a  given  time  frame.  The  time  frame  for  this  study  will  be 
one  day.  The  sortie  generation  rate  Is  a  function  of  the 
total  number. of  aircraft  at  a  location,  aircraft 
availability,  mission  length,  and  flight  line  operations. 

Statement  of  the  Problem 

The  problem  this  study  considers  Is  determining  the 
expected  aircraft  availability  and  sortie  generation  rate  for 
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the  T-46A  aircraft.  This  will  be  accomplished  by  simulating 
the  scheduled  and  unscheduled  maintenance  of  the  T-46A 
aircraft  as  well  as  the  flight  line  operations  at  a  typical 
Undergraduate  Pilot  Training  (OPT)  base. 

Scheduled  maintenance  Is  maintenance  that  Is  performed 
on  all  aircraft  after  a  specified  amount  of  flight  hours  have 
been  accumulated.  An  example  of  scheduled  maintenance  Is  the 
removal  and  thorough  Inspection  of  the  engines  after  every 
300  hours  of  flight  time.  Unscheduled  maintenance  is 
maintenance  that  Is  required  whenever  a  component  or 
subsystem,  such  as  a  radio,  falls  on  the  aircraft. 

The  flight  line  operations  Include  preflight  and 
postflight  Inspections,  on-aircraft  maintenance,  engine  test 
facilities,  engine  repair  shops,  and  other  off-alrcraft 
maintenance  activities.  This  study  will  be  limited  to 
preflight  and  postflight  Inspections  and  on-alrcraft 
maintenance.  Because  of  this  limitation,  this  study  will  not 
attempt  to  determine  manpower  requirements  at  the  UPT  base. 

The  model  will  be  based  upon  estimated  reliability  and 
maintainability  data  In  order  to  predict  the  expected 
performance  of  the  T-46A  aircraft  prior  to  actual  operational 
testing.  As  operational  testing  Is  begun  In  April  1986,  the 
test  data  that  Is  collected  will  be  used  to  update  the  model 
and  reevaluate  the  performance  characteristics  to  determine 
whether  the  T-46A  Is  achieving  the  desired  levels  of 
performance . 


Research  Question 


What  levels  of  aircraft  availability  and  sortie 
generation  rate  can  be  achieved  with  the  T-46A  trainer  In  an 
UPT  environment? 

Ob jec  tl ves 

The  overall  objective  of  this  research  is  to  develop  a 
valid  simulation  model  of  a  UPT  wing  equipped  with  T-46A 
aircraft  In  order  to  determine  the  expected  aircraft 
availability  and  sortie  generation  rates. 

In  order  to  fulfill  this  objective,  several 
subobjectives  need  to  be  accomplished.  These  subobjectives 
are  : 

1.  Collect  data  on  break  rates  and  repair  times  for 
major  subsystems  of  the  T-46A  aircraft. 

2.  Model  the  flying  operation  and  maintenance  of  the 
T-46A  aircraft  In  the  UPT  environment. 

3.  Develop  confidence  intervals  and  prediction 
equations  for  model  estimates  of  aircraft 
availability  and  sortie  generation  rates. 

General  Technique 

The  general  technique  that  will  be  used  In  this  research 
Is  simulation.  Simulation  Is  chosen  over  an  analytic 
technique  because  of  the  probabilistic  nature  of  modeling 
aircraft  flying  operations.  The  overall  reliability  and 
maintainability  of  an  aircraft  is  dependent  on  many  random 
processes.  These  random  processes  often  Interact  with  each 
other  which  makes  the  problem  of  determining  the  availability 


and  sortie  generation  rate  very  difficult  to  solve 


analytically.  Simplifying  assumptions  can  be  made  to  make 
the  problem  analytically  tractable,  however,  these  numerous 
assumptions  may  cast  doubt  on  the  validity  of  the  results, 
simulation,  on  the  other  hand,  can  model  the  interactions 
between  random  processes  and  provide  valid  results. 

The  simulation  will  be  accomplished  using  SLAM 
(Simulation  Language  for  Alternative  Modeling).  SLAM  Is 
chosen  for  several  reasons.  First,  the  sponsor  for  this 
research  (AFOTEC)  has  requested  that  SLAM  be  used.  AFOTEC 
will  be  using  this  simulation  model  during  operational 
testing  and  they  are  familiar  with  the  SLAM  language.  In 
addition,  SLAM  is  a  very  flexible  simulation  language  yet  it 
is  also  very  easy  to  use.  Finally,  a  SLAM  model,  being 
FORTRAN  based,  is  easily  transported  between  different 
computers.  Thus,  AFOTEC  will  be  able  to  easily  adapt  the 
model  to  their  computer. 

Methodology 

As  just  discussed,  simulation  is  the  general  technique 
chosen  for  providing  answers  to  the  research  question. 
However,  the  actual  methodology  required  to  arrive  at  those 
answers  Involves  accomplishing  the  subobjectives  that  were 
mentioned  earlier:  collecting  data,  modeling  the  T-A6A  in 
the  UPT  environment,  and  developing  confidence  Intervals  for 
the  model  output. 

Collecting  data  is  the  initial  phase  of  this  research. 
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Data  must  be  collected  which  estimates  the  break  rates  and 
repair  times  of  major  subsystems  of  the  T-46A  aircraft. 

Since  the  T-46A  is  a  new  aircraft,  past  performance  data  for 
the  subsystems  are  not  available.  However,  some  of  the 
subsystems  are  similar  to  subsystems  currently  being  used  on 
other  Air  Force  aircraft.  In  these  instances,  the  data  from 
the  currently  used  subsystems  can  be  substituted  for  the 
T-46A  subsystems.  The  reliability  for  the  subsystems  where 
comparability  data  does  not  exist  has  to  be  estimated. 
Knowledgeable  personnel  from  the  contractor  and  ATC 
maintenance  can  provide  realistic  estimates  for  these 
subsys  terns . 

After  collecting  the  data,  the  next  phase  is  to  model 
the  T-46A  aircraft  in  the  UPT  environment.  However,  before  a 
model  can  be  built,  the  modeler  must  understand  the  real  life 
system.  Therefore,  the  first  priority  is  gathering 
information  on  the  T-37  operations  in  the  current  UPT 
environment.  Once  the  current  system  is  understood,  the  next 
step  is  to  attempt  to  model  the  UPT  environment  as  accurately 
as  possible.  While  developing  the  model,  the  T-37  operations 
should  be  modified  to  reflect  anticipated  changes  in  the 
flight  operations  for  the  T-46A.  An  important  step  in  this 
phase  is  the  verification  and  validation  of  the  model  as  it 
is  built. 

Once  the  final  model  has  been  verified  and  validated, 
the  last  phase  of  the  research  is  to  develop  confidence 


Intervals  for  aircraft  availability  and  sortie  generation 
rate.  This  can  be  accomplished  by  using  regression  analysis 
techniques  which  estimate  the  effects  of  independent 
variables  on  the  response  variables.  Performing  the 
regression  analysis  requires  that  several  tasks  be 
accomplished.  The  first  task  is  to  identify  factors  in  the 
model  which  could  influence  the  main  performance 
characteristics,  aircraft  availability  and  sortie  generation 
rate.  The  second  task  is  to  design  an  experiment  to  collect 
output  data  from  the  model.  The  next  task  is  to  perform  a 
regression  analysis  on  the  data  to  estimate  the  parameters  of 
the  functional  relationship  suggested  by  the  data.  Finally, 
the  last  task  is  to  perform  analysis  on  the  factors  so  that 
confidence  intervals  can  be  constructed  around  the  model 
estimates  of  aircraft  availability  and  sortie  generation 
rate. 

Scope 

The  scope  of  this  study  has  been  limited  in  three  areas. 


First,  this  analysis  will  use  only  one  scenario.  This 
scenario  represents  the  current  UPT  operations  modified  to 
reflect  expected  changes  for  the  T-46A.  Second,  the  analysis 
is  limited  to  on-alrcraft  maintenance  only.  Finally,  the 
analysis  is  limited  to  UPT  operations  of  only  the  T-46A.  The 
analysis  will  not  include  any  T-38  considerations. 


Overview 


The  remainder  of  this  thesis  consists  of  four  chapters. 
Chapter  II  gives  a  verbal  description  of  the  model,  detailing 
the  UPT  environment  and  the  major  assumptions  used  In 
building  the  model. 

Chapter  III  analyzes  the  factors  and  outlines  the 
experimental  design.  In  addition,  the  steps  required  to 
ensure  valid  simulation  output  are  discussed. 

Chapter  IV  provides  the  results  of  the  experimental 
design  and  the  analysis  of  those  results. 

The  final  chapter.  Chapter  V,  presents  the  confidence 
Intervals  for  the  main  performance  characteristics,  aircraft 
availability  and  sortie  generation  rate,  as  well  as 
discussing  the  conclusions  reached  during  the  course  of  this 


II.  Model  Description 


A  full  description  of  the  model  of  the  T-46A  In  the  UPT 
environment  requires  that  a  description  of  the  aircraft  and 
the  UPT  environment  be  presented  first.  Because  the 
simulation  language  used  has  an  important  Impact  on  the 
development  of  the  model,  a  brief  description  of  SLAM  will 
also  be  presented  before  the  description  of  the  model.  These 
will  be  followed  by  an  overview  of  the  model  and  then  a  more 
detailed  narrative  description  of  the  model.  The  assumptions 
Inherent  In  the  model  will  be  presented  next.  Finally,  the 
chapter  will  conclude  with  a  description  of  the  steps  taken 
to  verify  and  validate  the  model. 

T-46 A  Description 

The  T-46A  Is  a  twin  engine  aircraft  with  side  by  side 
seating.  In  Its  role  as  the  primary  phase  trainer  for  UPT, 
the  T-46A,  like  the  T-37B,  must  be  capable  of  performing 
several  different  training  missions.  The  cockpit  of  the 
T-46A  will  be  pressurized  and  contain  more  modern  avionics. 
These  and  other  characteristics  will  Improve  the  flight 
training  capabilities  of  the  T-46A  as  compared  to  the  T-37B 
(2:3;  4:1).  The  availability  of  the  T-46A  will  depend 
heavily  on  the  reliability  of  Its  system  components.  This 
study  will  analyze  the  T-46A  by  classifying  the  aircraft  Into 
74  subsystems.  These  subsystems  are  listed  In  Table  B.l  of 
Appendix  B. 


The  T-46 A  and  the  UPT  System 

The  T-A6A  is  projected  to  be  assigned  to  six  ATC 
training  bases.  They  are  Columbus  AFB  MS,  Laughlin  AFB  TX, 
Reese  AFB  TX,  Vance  AFB  OK,  Williams  AFB  AZ ,  and  Mather  AFB 
CA  for  Undergraduate  Navigator  Training  (4:8).  Laughlin  AFB 
will  be  the  first  base  to  receive  the  T-46A  for  UPT.  Flight 
operations  are  essentially  identical  at  each  UPT  base.  Thes 
operations  are  described  next. 

The  flight  operations  of  the  T-46A  in  UPT  can  be 
described  under  two  broad  headings,  flying  and  maintenance. 
There  are  two  categories  of  maintenance,  scheduled 
maintenance,  which  are  preventative  actions,  and  unscheduled 
maintenance,  which  are  corrective  actions.  These  will  be 
discussed  later.  The  dally  flying  use  of  the  T-46A  includes 
four  activities.  They  are  scheduling  which  aircraft  will 
fly,  preparing  the  scheduled  aircraft  for  that  days  flying, 
the  flights  themselves,  and  after  each  flight,  a  short 
Inspection  and  servicing  of  the  aircraft.  These  four 
activities  are  discussed  in  more  detail  next. 

Scheduling  aircraft  to  fly  on  a  particular  day  is 
accomplished  during  the  previous  night.  Not  all  aircraft 
will  be  flown  each  day.  Some  will  not  be  available  to  fly 
because  of  maintenance  requirements.  Moreover,  in  general. 


there  are  more  aircraft  available  than  are  needed  to  fly  a 
days  training  schedule.  Aside  from  Insuring  that  there  are 
enough  aircraft  available  to  fly  the  days  missions,  the  main 


objective  of  determining  which  aircraft  to  use  is  to  provide 
an  even  flow  of  aircraft  into  the  phase  Inspection.  A  phase 
inspection  is  required  for  an  aircraft  after  it  has  been 
flown  a  specific  number  of  hours.  Therefore,  this  goal 
translates  to  keeping  the  number  of  flying  hours  all  aircraft 
have  been  flown  since  their  last  phase  inspection  uniformly 
distributed  from  zero  to  the  number  of  flying  hours  at  which 
a  phase  inspection  is  required  (14).  The  scheduling 
procedure  also  identifies  a  few  aircraft  to  be  used  as  spares 
in  case  the  primary  aircraft  are  unable  to  fly  during  the 
day.  After  the  scheduling  procedure  has  determined  which 
aircraft  will  be  flown  they  are  prepared  for  that  days 
missions. 

The  T-46A  will  receive  one  major  inspection  every  24 
hours  in  preparation  for  flying.  The  inspection  will  be  done 
in  place  of  the  two  Inspections,  one  before  the  days  flying 
and  one  after,  that  are  currently  being  done  for  the  T-37B 
(13).  This  daily  inspection,  hereafter  called  the  preflight 
inspection,  is  accomplished  during  the  early  morning.  The 
preflight  inspection  is  done  by  crew  chiefs  beginning  duty  at 
midnight  for  just  that  purpose.  These  crew  chiefs  will  also 
do  minor  maintenance  on  the  aircraft  if  needed. 

Once  the  aircraft  have  been  prepared,  they  may  fly 
several  times.  Most  flights  are  performed  during  daylight. 
These  flights  begin  15  minutes  before  sunrise  and  must  be 
finished  by  15  minutes  after  sunset.  If  required  for 


training,  night  flights  are  also  performed.  All  flights  are 
scheduled  with  a  minimum  of  three  minutes  between  takeoffs. 
Each  aircraft  on  the  schedule  may  be  flown  four,  five,  or  six 
times.  The  number  of  flights  an  Individual  aircraft  files 
depends  on  that  days  flying  schedule.  The  flying  schedule 
Includes  such  factors  as  the  amount  of  daylight  hours  and  the 
need  for  any  night  time  flying  training.  After  a  flight, 
each  aircraft  is  serviced.  The  servicing  includes  refueling 
and  a  short  walk-around  Inspection,  called  the  thrufllght 
Inspection.  In  addition  to  these  flying  activities,  the 
T-46A  also  undergoes  maintenance  actions.  These  actions  will 
be  described  next. 

Maintenance  actions  required  by  the  T-46A  fall  Into  two 
categories,  scheduled  and  unscheduled.  Scheduled  maintenance 
Is  preventative  maintenance  and  is  to  be  done  to  keep  the 
aircraft  in  a  ready- to-f ly  status.  This  maintenance  Includes 
the  preflight  and  thrufllght  inspections,  corrosion 
prevention,  and  phase  inspections,  among  other  scheduled 
maintenance  actions.  Scheduled  maintenance  Is  required  based 
on  the  number  of  hours  the  aircraft  has  been  flown.  Table 
B.3  of  Appendix  B  contains  a  list  of  scheduled  maintenance 
actions  and  the  flying  hour  Intervals  between  them  as 
proposed  by  Fairchild  Republic.  Scheduled  maintenance  Is 
performed  mainly  during  the  normal  dayshift. 

Unscheduled  maintenance  Is  corrective  maintenance  done 
to  return  an  aircraft  to  a  ready-to-fly  status  after  a  part 


has  failed  or  has  been  reported  as  malfunctioning. 

Unscheduled  maintenance  Is  performed  when  needed  but  the 
majority  of  it  Is  performed  during  the  swing  shift,  from  1600 
to  midnight  (14).  These  maintenance  actions  are  performed  by 
technicians  assigned  by  specialties  to  twenty  different  work 
centers.  Table  B.2  of  Appendix  B  contains  a  list  of  the  work 
centers  currently  used  in  the  ATC  maintenance  policy. 

In  summary,  there  are  two  major  activities  in  the 
operation  of  the  T-46A  in  its  role  as  the  primary  trainer 
aircraft  in  UPT.  They  are  flying  activities  and  maintenance 
activities.  The  flying  activities  include  determining  which 
aircraft  will  fly  on  a  particular  day,  preparing  them  for 
flight,  the  flying  itself,  and  post  flight  servicing.  The 
maintenance  activities  include  scheduled  and  unscheduled 
maintenance.  Before  describing  the  model  of  this  system,  a 
background  description  of  SLAM  is  necessary  to  understand  how 
this  model  was  developed. 

SLAM  Background 

SLAM  is  a  special  purpose  language  which  is  used  for 
simulation  modeling.  It  is  based  on  the  FORTRAN  language. 
SLAM  provides  two  orientations,  or  a  combination  of  both, 
to  modeling.  They  are  e ven t- schedu 1 1 ng  and 
process-interaction  (6:99).  Each  orientation  has  its 


advantages  and  disadvantages.  The  pr oce s s- in te r ac 1 1  on 
orientation  is  easier  to  use  but  may  not  describe  all  the 
processes  that  can  occur.  The  event-scheduling  orientation 


allows  modeling  to  the  level  of  complexity  desired  but  at  a 
cost  of  an  increased  modeling  effort  (12:323).  Fortunately, 
SLAM  allows  both  orientations  to  be  used  simultaneously. 

The  process-interaction  orientation  of  SLAM  uses 
networking  concepts  to  model  a  system.  There  are  nodes  and 
branches  which  represent  parts  of  the  system  such  as  decision 
points,  queues,  and  activities.  Entities,  such  as  aircraft 
in  this  case,  then  flow  through  the  network. 

The  event-scheduling  orientation  of  SLAM  uses  the 
concept  that  changes  to  the  system  can  be  modeled  as 
happening  at  specific  Instances.  These  changes,  called 
events,  are  coded  in  FORTRAN  subroutines  by  the  modeler. 

These  events  can  be  as  complex  as  needed  to  model  the  system. 
SLAM  automatically  controls  time  advancement  and  sequencing 
of  events.  SLAM  also  provides  subprograms  that  can  be  used 
for  common  event  activities  such  as  event  scheduling,  random 
sampling,  and  statistics  collection  (12:73). 

The  SLAM  model  developed  for  the  T-46A  uses  both  the 
pr oce s s- 1 n te r ac t i on  and  e ve n t-schedu 1 i ng  orientations.  Use 
of  both  orientations  allows  for  entitles  in  the  network  model 
to  initiate  events  and  for  events  to  change  the  flow  of 
entitles  in  the  network  (12:74).  An  overview  of  the  model  is 
provided  before  a  more  detailed  narrative  description  of  the 
model  is  presented. 

Model  Overview 

Earlier,  the  UPT  system  was  described  as  consisting  of 


two  major  activities,  flying  activities  and  maintenance 
activities.  The  purpose  of  the  model  overview  is  to  present 
the  model  structure  and  show  how  the  flying  and  maintenance 
activities  are  incorporated  into  this  structure.  The 
description  of  the  model  structure  will  include  an 
explanation  of  why  the  repair  network  is  not  a  SLAM  network 
and  a  description  of  how  the  aircraft  ate  modeled.  In 
addition  to  the  model  structure,  the  overview  will  present 
information  about  the  data  for  the  model.  This  will  Include 
the  number  of  aircraft  chosen,  manpower  resource  levels  used, 
and  the  data  sources  for  the  model. 

Model  Structure .  The  T-46A  model  is  a  combined  network 
and  discrete-event  simulation  model.  It  consists  of  two 
parts,  a  SLAM  network  portion  and  a  FORTRAN  portion.  The 
SLAM  portion  consists  of  three  major  network  segments  and 
four  supporting  network  modules.  The  three  major  network 
segments  are  the  sortie  generation,  failure,  and  phase 
inspection  segments.  The  sortie  generation  segment  includes 
all  four  actions  described  as  flying  activities.  The  failure 
segment  covers  unscheduled  maintenance  while  the  phase 
inspection  segment  provides  for  scheduled  maintenance.  These 
three  major  network  segments  are  Interconnected.  The  sortie 
generation  segment  includes  branches  to  the  other  two  major 
segments  at  the  appropriate  times.  Figure  1  shows  the 
relationship  of  the  major  model  segments. 

The  interconnection  of  the  major  network  segments  is  in 
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contrast  to  the  four  supporting  network  modules.  Each  of 
these  is  an  Independent  network.  However,  the  four  support 
modules  control  the  flow  of  activities  in  the  three  major 
network  segments.  These  modules  limit  flying  to  daylight 
only,  change  crew  sizes  at  shift  changes,  create  three  minute 
intervals  between  takeoffs,  and  reset  counters  used  for 
statistics. 

The  remaining  portion  of  the  model  is  the  FORTRAN 
program.  The  FORTRAN  program  consists  of  two  major  parts,  an 
allocation  subroutine  and  an  event  subroutine.  The 
allocation  subroutine  allocates  repair  crew  resources  to  fix 
aircraft  with  failures.  The  event  subroutine  controls  the 
discrete-event  orientation  of  the  model.  One  part  of  the 
event  subroutine  is  the  repair  network. 

The  repair  network  for  unscheduled  maintenance  is  based 
on  the  ATC  Logistics  Composite  Modeling  (LOOM)  model  network. 
However,  because  of  a  SLAM  limitation  of  300  nodes,  as 
Implemented  on  the  AFIT  VAX  11/785  computer,  the  LOOM 
unscheduled  maintenance  network  cannot  be  converted  to  a  SLAM 
network.  Thus,  the  network  is  contained  in  a  FORTRAN 
subroutine  and  is  somewhat  hidden  from  view.  The  form  of 
this  network  will  be  discussed  in  more  detail  in  the 
narrative  description  of  the  model.  There  are  two  additional 
reasons  for  choosing  a  FORTRAN  network.  The  first  is  that 
entering  or  changing  data  for  the  FORTRAN  network  is  easier 
than  for  a  SLAM  network.  The  second  and  more  Important 


reason  is  that  FORTRAN  allows  easier  modification  of  the 
repair  network.  For  example,  consider  a  change  in  the  work 
unit  code  (WUC)  level  at  which  the  aircraft  is  modeled.  The 
WUC  is  an  Indication  of  the  amount  of  detail  at  which  the 
aircraft  is  analyzed.  The  number  of  digits  in  the  WUC 
corresponds  to  the  level  of  aggregation.  For  example,  a 
5-digit  WUC  represents  a  specific  part  whereas  a  3-dlglt  WUC 
represents  a  subsystem  such  as  the  nose  landing  gear.  For  an 
Increase  in  the  number  of  WUC  digits,  the  FORTRAN  network 
requires  only  a  change  in  the  global  variable  defining  the 
number  of  WUCs  used.  SLAM,  on  the  other  hand,  requires  that 
a  new  network  segment  be  built  for  each  new  WUC.  The  level 
of  detail  chosen  for  this  model  is  presented  next. 

In  determining  the  overall  availability  of  the  T-46A  the 
aircraft  was  modeled  at  the  3-digit  WUC  level.  There  are  74 
3-digit  WUCs  modeled.  Table  B.l  of  Appendix  B  contains  a 
list  of  the  WUCs  and  the  mean  time  between  failures  (MTBF) 
for  each.  The  failures  of  these  subsystems  are  generated 
using  a  probability  per  flight  rather  than  a  failure  clock. 

So  far,  the  overview  has  presented  the  structure  of  the 
model.  The  remainder  of  the  overview  will  present 
information  about  the  data  used  in  the  model. 

Data  Input.  Since  Laughlin  AFB  is  scheduled  to  receive 
the  T-46A  first,  the  number  of  aircraft  and  manpower  levels 
from  that  base  will  be  used  as  representative  levels. 

Laughlin  currently  has  82  T-37B  aircraft  assigned  to  it  for 


UPT.  This  is  the  number  of  T-46A  aircraft  used  In  modeling 


this  UPT  system.  The  manpower  levels  of  repair  technicians 
can  be  found  In  Table  B.2  of  Appendix  B.  It  Is  presented 
there  by  work  center  and  shift. 

There  are  three  major  sources  of  data  for  this  model. 
Information  concerning  the  T-46A  in  the  UPT  system  came  from 
Interviews  conducted  with  MAJ  Schad  (14),  Quality  Assurance 
at  Laughlln  AFB,  and  MAJ  Purcell  (13),  T-46A  System  Program 
Office  at  Wright-Patterson  AFB,  OH.  The  ATC  LOOM  Final 
Report  (3)  and  model  (5)  provided  numerous  data  on  manpower 
levels,  the  UPT  system,  and  the  unscheduled  maintenance 
network.  An  early  version  of  the  Aeronautical  System 
Division  (ASD)  T-46A  LOOM  model  (1)  supplemented  the  ATC  LOOM 
model  on  the  last  Item.  Finally,  a  Fairchild  Republic 
document  (8)  provided  estimates  of  the  MTTR  (mean  time  to 
repair)  and  MTBF  of  the  subsystems  of  the  T-46A.  This  report 
also  provided  Information  about  scheduled  maintenance.  The 
data  can  be  found  In  Appendix  B. 

Narrative  Description 

The  FORTRAN  portion  of  the  model  will  be  described  first 
because  the  FORTRAN  code  defines  events  that  are  used  In  the 
SLAM  network.  The  SLAM  and  FORTRAN  codes  are  contained  In 
Appendix  A. 

The  FORTRAN  model  consists  of  five  parts,  the  main 
program  which  Initiates  the  simulation,  an  Initialization 
subroutine,  an  event  subroutine,  an  allocation  subroutine. 
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and  an  output  subroutine.  Only  the  event  and  allocation 
subroutines  will  be  described  here. 

There  are  nine  events  that  are  accessed  through  the 
event  subroutine.  They  are  named  DAYSHIFT,  SWINGSHIFT, 
MIDSHIFT,  SCHEDULE,  TOPREFLIGHT,  TAKEOFF,  FAIL,  FREECREWS  , 
and  DETCREWS.  Three  of  the  events  model  the  change  In  the 
number  of  repair  personnel  at  shift  changes.  These  are  the 
events  DAYSHIFT,  SWINGSHIFT,  and  MIDSHIFT.  Event  SCHEDULE 
determines  which  aircraft  will  be  flown  the  next  day.  This 
event  also  causes  event  TOPREFLIGHT  to  occur.  This  Is  the 
event  which  actually  allows  the  aircraft  chosen  by  SCHEDULE 
to  be  queued  for  a  prefllght.  The  event  TAKEOFF  allows  one 
aircraft  which  Is  ready  to  fly  to  begin  the  flight  sequence. 
Event  FAIL  determines  If  any  of  the  aircrafts  subsystems  have 
failed.  This  Is  done  by  converting  the  MTBF  of  each 
subsystem  to  a  probability  of  failure  for  one  flight.  Then  a 
random  number  from  zero  up  to  one  Is  drawn  and  checked  to  see 
If  It  falls  In  the  failure  range  for  that  subsystem.  All  of 
the  aircraft  systems  are  checked  for  a  failure.  The  event 
FREECREWS  releases  crews  which  have  been  used  to  repair  an 
aircraft.  The  final  event  Is  DETCREWS  which  Is  the 
unscheduled  repair  network.  Event  DETCREWS  determines  which 
crew  will  repair  a  failure.  The  size  of  the  repair  crew  Is 
also  determined.  If  an  aircraft  has  more  than  one  failure, 
the  crews  needed  to  repair  all  failures  are  determined  at 
this  time.  Event  DETCREWS  models  the  repair  network  in  a  way 


similar  to  LCOM.  The  ATC  LCOM  model  of  the  T-37B  provided 
the  basic  form  for  this  model's  repair  network.  However,  the 
T-46A  has  a  subsystem  structure  which  is  slightly  different 
than  the  T-37B.  Therefore,  an  early  version  of  the  ASD  LCOM 
model  of  the  T-46A  was  used  to  supplement  the  ATC  LCOM  model. 
This  T-46A  repair  network  model  Includes  repair  activities 
that  have  one  or  two  actions.  The  term  action  is  used  here 
to  represent  the  steps  necessary  to  complete  an  LCOM  repair 
task.  In  LCOM,  tasks  are  coded  into  several  categories  such 
as  minor  maintenance  ("M")  or  cannot  duplicate  ("H"). 

However,  for  the  tasks  taken  from  the  ATC  network  model,  the 
"M"  tasks  consolidate  minor  maintenance  and  cannot  duplicate 
tasks.  The  "R"  tasks  denote  remove  and  replace  tasks  (3:1-1 
to  1-2). 

Most  of  the  repair  activities  in  this  model  have  only 
one  maintenance  action  and  that  action  may  be  performed  by 
one  of  a  few  specific  repair  crews.  A  few  maintenance 
activities  have  two  actions.  The  first  action  is  one  that 
must  be  done  first  by  a  specific  crew  but  the  other  may  be 
performed  by  one  of  several  crews.  An  example  of  a 
maintenance  activity  with  two  actions  is  repairing  a  main 
landing  gear.  The  aircraft  must  first  be  put  on  jacks  before 
the  actual  repair  activity  takes  place.  In  this  model  those 
actions  which  must  be  done  prior  to  the  actual  repair 
activity  are  called  required  maintenance  actions.  Required 
maintenance  actions  are  always  done  by  the  one  type  of  repair 
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crew.  All  other  maintenance  actions  are  called  possible 
maintenance  actions  because  one  of  several  different  types  of 
crews,  possibly  with  different  si2es,  may  complete  the  repair 
activity.  The  repair  network  is  contained  in  Table  B.4  of 
Appendix  B.  The  other  major  portion  of  the  FORTRAN  program 
is  the  allocation  subroutine. 

The  allocation  subroutine  seizes  the  repair  crew  needed 
to  fix  an  aircraft  when  the  repair  crew  is  available.  If  an 
aircraft  has  more  than  one  failure,  the  availability  of  all 
necessary  repair  crews  are  checked.  If  there  is  more  than 
one  aircraft  waiting  for  repair  crews,  all  of  the  aircraft 
are  checked  to  see  if  crews  are  available  to  fix  them.  A 
more  detailed  description  of  the  allocation  subroutine  can  be 
found  in  Appendix  A. 

Figure  2  shows  the  sortie  generation  network  in  detail. 
First,  32  aircraft  are  created  and  assigned  an  initial  number 
of  flight  hours.  The  initial  flight  hours  are  distributed 
uniformly  from  0  to  300.  Three  hundred  hours  is  the  interval 
Fairchild  Republic  recommends  between  phase  inspections 
(3:4-10).  The  aircraft  then  wait  to  be  scheduled  to  fly  the 
next  days  missions.  Approximately  40  aircraft  are  normally 
used  on  one  day  (14).  However,  in  this  model  all  aircraft 
are  allowed  to  fly  in  order  to  determine  what  sortie 
generation  rate  the  T-46A  can  achieve.  The  aircraft  which 
have  been  chosen  to  fly  are  then  given  a  preflight 
inspection.  The  aircraft  next  wait  for  daylight  and  a 
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takeoff  slot.  The  length  of  the  flight  is  randomly  set  based 
on  a  normal  distribution  with  a  mean  of  1.3  hours  and  a 
variance  of  10%  (3:5-3).  Following  the  flight,  the  aircraft 
Is  sent  to  the  phase  Inspection  network  if  the  aircraft  has 
accumulated  more  than  300  hours.  When  a  phase  inspection  is 
not  required  the  aircraft  is  then  sent  to  the  failure 
network.  The  aircraft  may  return  from  the  failure  network  to 
two  different  places  in  the  sortie  generation  network.  If 
more  than  twenty-four  hours  have  passed  since  its  last 
preflight,  the  aircraft  is  sent  to  be  scheduled  for 
the  next  day.  When  a  preflight  is  not  required,  the  aircraft 
is  sent  back  to  a  decision  node  called  CONT  in  the  sortie 
generation  network.  From  here  the  aircraft  is  sent  to  wait 
for  a  takeoff  slot  if  it  is  still  daylight.  If  it  is  not, 
the  aircraft  is  sent  to  be  scheduled  for  the  next  day.  At 
the  end  of  the  day,  all  aircraft  that  are  waiting  for  takeoff 
are  also  sent  to  be  scheduled  for  the  next  day. 

When  a  phase  Inspection  is  required,  the  aircraft  first 
waits  for  a  phase  dock  to  become  available.  There  are  four 
phase  docks  in  this  model.  Maintenance  manhours  for  all 
scheduled  maintenance  actions  are  accumulated  in  this  network 
with  the  exception  of  those  hours  spent  on  preflight  and 
thruflight  inspections.  These  maintenance  manhours  are 
collected  as  they  occur  in  the  sortie  generation  network. 
Following  completion  of  the  phase  inspection  the  flying  hours 
for  that  aircraft  are  reset  to  zero  and  the  aircraft  is  sent 
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to  be  scheduled  for  the  next  days  missions. 


Figure  3  contains  a  diagram  of  the  failure  network. 
First,  each  aircraft  is  checked  to  determine  if  a  failure  has 
occurred.  If  there  are  no  failures  the  aircraft  is  sent  back 
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Figure  3.  Failure  Network  Segment 


to  the  sortie  generation  network.  When  there  are  failures, 
the  crews  needed  to  repair  the  failures  are  determined.  The 
aircraft  then  waits  for  these  crews  to  become  available  to 
perform  the  unscheduled  maintenance.  After  all  failures  have 


been  fixed,  the  aircraft  is  returned  to  the  sortie  generation 
network.  The  aircraft  may  return  to  that  network  In  two 
different  places.  If  the  aircraft  has  had  a  prefllght  within 
the  last  24  hours  It  Is  returned  where  It  can  be  flown  If 
there  Is  still  daylight.  When  an  aircraft  needs  a  preflight 
It  Is  sent  to  be  scheduled  for  the  next  day. 

This  T-46A  model  has  been  developed  to  a  degree 
sufficient  to  test  the  effects  of  various  factors  on  the 
availability  of  the  T-46A  and  its  sortie  generation  rate. 
There  are,  however,  assumptions  Inherent  In  the  model  that 
effect  the  prediction  of  availability  and  sortie  generation 
rate.  These  are  discussed  next. 

Model  Assumptions 

The  assumptions  Inherent  In  this  model  fall  Into  two 
categories.  Assumptions  created  by  leaving  something  out  of 
the  model  and  assumptions  made  In  determining  the  working 
details  of  the  model.  In  the  former  category,  this  model 
does  not  address  spare  parts  or  weather.  The  Inclusion  of 
spare  parts  Is  beyond  the  scope  of  this  study.  Resupply  and 
cannibalization  of  spares  are  complex  Issues  which  may 
warrant  further  study  but  they  are  not  considered  here.  A 
crude  approximation  of  the  effect  of  spare  parts  on 
availability  can  be  made  by  subtracting  the  historic 
percentage  of  aircraft  not  mission  capable  due  to  supply. 

This  reduction  In  the  number  of  aircraft  available  is 
anticipated  to  have  no  effect  on  the  T-46A  sortie  generation 
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rate.  Weather  was  not  modeled  because  the  Intent  of  this 
study  Is  to  determine  what  the  sortie  generation  rate  of  the 
T-46A  will  be  If  allowed  to  fly  under  the  ATC  three  minute 
stratified  takeoff  concept.  ATC  considers  the  Impact  weather 
will  have  on  the  UPT  system  when  determining  the  flying 
training  programs  for  each  base.  Thus,  there  Is  more 
Interest  on  what  the  aircraft  can  do  when  not  constrained  by 
the  weather. 

There  are  eight  assumptions  made  In  determining  the 
logic  of  the  model.  They  are: 

1.  A  repair  crew  works  until  finished  with  a  repair. 

2.  Most  scheduled  maintenance  manhours  are  counted  In  the 
phase  Inspection  network. 

3.  Scheduling  of  aircraft  takes  place  at  midnight. 

4.  An  aircraft  Is  checked  for  failures  after  each  flight. 

5.  Multiple  failures  are  repaired  sequentially. 

6.  Multiple  failures  are  repaired  from  lowest  WUC  to 
hlghes  t. 

7.  Aircraft  waiting  for  a  particular  crew  are  repaired  In 
random  order. 

8.  Aircraft  are  given  prefllght  Inspections  In  order  of 
flight  hours. 

The  first  assumption  made  In  the  logic  of  the  model  Is 
that  once  a  repair  Is  started  the  repair  crew  will  work  until 
It  Is  finished  with  the  repair.  This  assumption  may  have  an 
Impact  when  the  crew's  shift  Is  scheduled  to  end  during  the 
repair.  The  Impact  Is  felt  only  when  the  shift  change 
results  In  a  decrease  of  repair  personnel  on  duty  for  that 


specialty  code.  In  this  case,  the  decision  to  keep  the  crew 
working  would  normally  be  made  by  the  supervisor.  The 
decision  would  be  based  on  the  need  for  that  particular 
aircraft  and  other  considerations.  However,  for  repairs  that 
last  just  a  little  past  a  shift  change  the  impact  of  this 
assumption  is  negligible. 

A  second  assumption  involves  how  the  manhours  for 
scheduled  maintenance  are  counted.  Scheduled  maintenance 
Includes  those  actions  listed  in  Table  B.3  and  also  the 
preflight  and  thrufllght  inspections.  The  maintenance 
manhours  for  the  preflight  and  thrufllght  Inspections  are 
counted  as  these  actions  occur  in  the  sortie  generation 
network.  There  is  only  one  network  for  the  remaining 
scheduled  maintenance,  the  phase  inspection  network.  The 
manhours  for  all  maintenance  listed  in  Table  B.3  are  counted 
during  this  network.  This  is  done  by  calculating  a 
per-300-hour  equivalent  for  all  of  the  maintenance  in  Table 
B.3.  This  sura,  which  is  97.798,  is  then  counted  each  time  an 
aircraft  enters  phase.  It  may  be  more  realistic  to 
accumulate  more  of  the  maintenance  manhours  as  each  scheduled 
maintenance  action  should  occur. 

When  the  scheduling  activity  takes  place  is  another 
assumption.  This  is  a  concern  because  the  preflight 
inspections  can  begin  only  after  the  scheduling  is  complete. 
In  the  model,  scheduling  is  done  at  midnight  every  night. 

This  allows  the  midnight  shift  crew  chiefs  eight  hours  to 


preflight  aircraft.  In  the  real  world  scheduling  may  be 
finished  before  midnight  thus  allowing  more  time  to  preflight 
aircraft.  However,  eight  hours  for  the  midnight  shift  plus 
the  additional  time  in  the  morning  before  the  first  aircraft 
returns  from  flying  should  allow  all  aircraft  to  be 
preflighted  before  the  crew  chiefs  are  also  needed  for  repair 
activities. 

The  premise  that  an  aircraft  is  checked  for  failures 
after  each  flight  is  a  fourth  assumption.  This  does  not 
happen  when  the  aircraft  is  required  to  have  a  phase 
inspection.  In  this  case  the  aircraft  is  sent  to  the 
inspection  In  lieu  of  being  checked  for  failures.  Not 
checking  for  failures  after  every  flight  may  have  an  Impact 
on  the  model  because  the  failures  are  being  generated  based 
on  a  probability  of  failure  per  flight.  However,  the  impact 
is  expected  to  be  small.  If  the  impact  is  large,  it  will  be 
seen  in  a  chi-square  statistic  computed  to  test  the  validity 
of  the  failure  generator. 

Another  model  assumption  concerning  failures  is  that 
they  are  always  repaired  sequentially.  Some  maintenance 
activities  such  as  those  involving  repair  of  the  fuel  system 
are  not  allowed  to  be  done  concurrently  with  any  other 
repair.  Other  combinations  of  repairs  may  be  impractical. 
Since  data  on  which  repairs  are  allowed  concurrently  is  not 
available,  the  model  performs  repairs  sequentially.  Further, 
during  unscheduled  maintenance,  sequential  repairs  happen 


more  frequently  than  concurrent  repairs. 

The  remaining  three  assumptions  in  the  logic  of  the 
model  deal  with  the  order  in  which  aircraft  entitles  are 
moved  within  the  SLAM  network.  This  concern  about  the  order 
arises  twice  In  the  allocation  subroutine.  The  first  of 
these  probably  has  little  or  no  effect.  It  is  the  order  in 
which  multiple  failures  to  an  aircraft  are  repaired.  The 
model  checks  for  crew  availability  in  the  order  of  the  lowest 
numbered  failure  to  the  highest.  Because  of  the  way  failures 
are  determined,  the  failures  for  an  aircraft  are  ordered  from 
lowest  U(JC  to  highest.  The  first  failure  to  have  crews 
available  for  the  repair  is  fixed. 

The  second  concern  may  have  more  of  an  impact.  This 
concern  results  from  switching  the  order  of  the  aircraft 
waiting  for  repair  crews  when  there  Is  more  than  one  aircraft 
waiting.  The  switching  is  done  to  Insure  that  when  a  crew 
becomes  available  to  repair  a  waiting  aircraft  that  the 
repair  begins  at  that  time.  The  switching  allows  for  the 
possibility  that  when  two  aircraft  are  waiting  for  the  same 
crew,  the  aircraft  that  has  been  waiting  the  least  time  may 
be  repaired  first.  This  may  be  done  in  actual  practice  and, 
in  any  event,  should  not  effect  such  statistics  as  the 
average  waiting  time.  However,  it  will  effect  the  longest 
waiting  time . 

The  order  in  which  the  aircraft  are  placed  in  the  queue 
for  preflights  is  the  final  model  assumption.  In  the  model 


aircraft  are  placed  in  this  queue  during  the  scheduling 


event.  This  event  orders  the  aircraft  by  the  number  of 
flight  hours  and  then  places  them  In  the  preflight  queue  by 
that  ranking.  The  reason  for  the  concern  here  Is  that 
aircraft  which  are  preflighted  first  and  thus,  fly  first, 
have  a  higher  probability  of  flying  more  missions  In  a  given 
day  than  aircraft  that  are  Inspected  later.  Whether 
ascending  or  descending  order  Is  used  to  rank  the  aircraft 
can  have  a  significant  Impact.  If  aircraft  are  placed  Into 
the  queue  In  ascending  order,  the  distribution  of  flight 
hours  becomes  less  uniform  as  the  aircraft  with  less  hours 
fly  more  times  per  day  than  those  with  more  hours.  This 
eventually  causes  many  aircraft  to  reach  phase  Inspection  at 
nearly  the  same  time.  This  model  places  the  aircraft  In  the 
queue  from  highest  to  lowest  number  of  flight  hours  allowing 
the  distribution  of  flight  hours  to  remain  more  uniform. 

This  ordering  concern  may  not  be  as  much  of  a  problem  If  the 
utilization  rate  of  the  T-A6A  Is  constrained  or  less  than  all 
aircraft  are  allowed  to  fly. 

Verification  and  Validation  of  the  Model 

Verification  Is  the  process  of  determining  the  model 
works  as  Intended  while  validation  Is  the  process  of 
determining  the  model  accurately  portrays  the  real  system 
being  modeled  (12:10).  This  model  was  verified  through  the 
use  of  two  techniques,  trace  listings  and  summary  reports. 


Trace  listings  showing  how  the  aircraft  entitles  moved 


through  the  SLAM  networks  and  the  FORTRAN  subroutines  were 
generated.  The  trace  listings  revealed  that  the  aircraft 
moved  through  the  SLAM  network  as  intended.  The  following 
are  examples  of  network  flows  that  were  observed.  Aircraft 
were  given  a  preflight,  flew  several  times,  stopped  flying  at 
night,  and  then  began  this  cycle  again.  Phase  inspections 
were  completed  at  the  appropriate  times.  Aircraft  were 
checked  for  failures  and  routed  correctly  based  on  whether  a 
failure  had  occurred.  In  addition,  aircraft  with  multiple 
failures  were  sent  through  the  repair  cycle  until  all 
failures  were  fixed.  The  trace  listings  also  showed  that  the 
actions  in  the  discrete  events  were  occurring  properly. 
Examples  of  these  are  the  following.  Failure  and  crew 
determinations  occurred  correctly  for  the  random  numbers 
drawn.  Crews  were  allocated  correctly  based  on  the  number 
needed  and  the  number  available.  If  there  were  more  than  one 
aircraft  waiting  for  crews,  all  aircraft  were  checked  when  a 
crew  became  available.  Shift  changes  occurred  correctly  at 
the  appropriate  times.  Also,  the  trace  listings  showed  that 
housekeeping  details,  such  as  which  planes  have  what  failures 
and  which  crews  were  being  waited  for,  were  kept  correctly. 

The  SLAM  summary  reports  were  examined  for  indications 
of  problems  such  as  unexpected  queue  lengths  and  destruction 
or  creation  of  aircraft  entities.  In  particular,  the  queues 
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The  queue  for  the  thrufllght  always  remained  at  a  reasonable 
length.  However,  at  first,  the  queue  length  for  the  phase 
docks  was  excessive  because  of  the  way  the  aircraft  were 
ordered  during  the  scheduling  event.  When  the  ordering  of 
the  aircraft  was  switched  to  those  with  the  highest  number  of 
flight  hours  received  prefllghts  first,  the  queue  length  for 
the  phase  docks  was  negligible. 

The  SLAM  summary  reports  were  used  to  show  that  entitles 
at  branches  In  the  network  did  not  get  lost  nor  were  any 
extra  entitles  created.  For  example,  the  number  of  aircraft 
leaving  the  two  exit  points  of  the  failure  network  equaled 
the  number  entering  this  network.  Each  branch  node  was 
tested  to  make  sure  that  the  number  of  entitles  leaving  the 
node  was  equal  to  the  number  entering  the  node. 

One  additional  step  was  taken  to  verify  the  model.  A 
chi-square  statistic  was  computed  to  test  whether  the 
failures  were  being  distributed  across  all  WUCs  correctly. 

The  chi-square  statistic  for  a  test  run  of  48,000  hours  was 
90.922.  The  critical  value  for  a  =  0.01  and  73  degrees  of 
freedom  Is  104.01  (9:437).  Therefore,  the  hypothesis  that 
the  failure  generator  works  correctly  cannot  be  rejected. 

This  chi-square  statistic  was  also  computed  for  all  of  the 
experimental  runs.  The  value  of  the  statistic  ranged  from 
60.775  to  89.395  for  these  runs. 

Validation  of  the  model  Is  a  more  difficult  task. 
Ideally,  a  model  can  be  validated  by  using  historic  Inputs 


and  then  comparing  the  model  outputs  to  the  historic  outputs. 
Since  the  T-46A  is  a  new  aircraft  there  is  no  historic  data 
that  can  be  used.  Therefore,  our  validation  efforts  were 
aimed  at  considering  the  reasonableness  of  the  model  outputs 
to  the  given  Inputs.  The  observed  output  values  were 
determined  to  be  near  the  expected  values  for  such  measures 
as  the  total  number  of  failures  and  MTTR.  In  addition, 
changes  in  the  output  measures  occurred  in  ways  expected  as 
the  inputs  were  varied.  For  example,  the  availability  of  the 
aircraft  decreased  when  the  failure  rate  was  increased. 


III.  Methodology 

The  purpose  of  this  chapter  is  to  select  factors  to  be 
considered  in  an  experimental  design  and  then  to  select  the 
most  appropriate  experimental  design.  In  addition,  the  steps 
taken  to  ensure  valid  simulation  output  for  the  experimental 
design  are  discussed. 

Factor  Selection 

A  critical  step  in  deciding  which  experimental  design  to 
use  is  determining  which  factors  need  to  be  examined. 
Initially,  four  factors  were  considered:  phase  inspections, 
manpower,  mean  time  between  failures  (MTBF),  and  mean  time  to 
repair  (MTTR). 

In  looking  at  the  effects  of  phase  inspections, 
increasing  the  frequency  of  phase  in  the  model  would  only 
decrease  the  aircraft  availability  because  the  aircraft  is 
tied  up  in  phase  more  often.  This  decrease  in  availability 
is  in  contrast  to  the  real  world  in  which  an  increase  in 
availability  is  possible  if  the  increased  frequency  of  phase 
makes  the  aircraft  more  reliable.  Because  this  relationship 
could  not  be  quantified  it  was  decided  not  to  include  phase 
Inspection  as  a  factor. 

There  are  two  reasons  manpower  was  not  used  as  a  factor 
in  this  analysis.  The  first  reason  is  the  structure  of  the 
manpower  data.  The  manpower  data,  as  previously  mentioned. 


was  obtained  from  the  current  LCOM  model  for  the  T-37 


operations  at  Laughlln  AFB.  In  some  cases  however,  the  data 
represented  personnel  for  both  the  T-37  as  well  as  the  T-38 
operations.  In  these  cases.  It  was  not  possible  to  identify 
the  number  of  people  who  worked  on  the  T-38.  In  addition, 
there  are  many  ways  to  apportion  manpower  to  provide  a 
constant  availability  or  sortie  generation  rate.  Because 
there  is  no  unique  solution  and  also  no  clear  guidance  on  how 
to  reduce  manpower,  the  Laughlln  manpower  data  could  no  be 
reduced  to  reflect  only  T-46A  maintenance.  Therefore,  this 
data  would  tend  to  overstate  the  T-46A  manpower. 

The  second  reason  is  the  intent  of  this  analysis.  While 
it  may  be  possible  to  reduce  the  given  levels  of  manpower  in 
the  model,  this  model  was  not  designed  to  estimate  what  the 
actual  manpower  requirements  should  be. 

For  these  reasons,  it  was  determined  that  varying 
manpower  would  not  provide  a  basis  for  meaningful  analysis 
and  therefore  should  not  be  included  as  a  factor.  However, 
it  would  be  of  Interest  to  determine  the  maximum  aircraft 
availability  and  sortie  generation  rate  that  could  be 
achieved  if  manpower  is  not  constrained  at  all. 

One  of  the  current  Issues  in  the  Air  Force  is  the 
acquisition  of  systems  that  are  reliable  and  easily 
maintained.  As  was  stated  in  the  Introduction,  one  of  the 
objectives  for  the  T-46A  is  operational  and  support  cost 
savings  through  the  use  of  modern  technology.  This 


translates  directly  to  the  reliability  and  maintainability  of 


an  aircraft.  Therefore,  it  was  decided  that  the  experimental 
design  should  focus  on  the  reliability  and  maintainability 
parameters  (MTBF  and  MTTR)  of  the  T-46A  and  their  effects  on 
aircraft  availability  and  sortie  generation  rate. 


The  Design 

Regression  analysis  was  the  approach  chosen  to 
investigate  the  effects  of  the  two  independent  variables, 

MTBF  and  MTTR,  on  the  response  variables.  Regression 
analysis  combines  an  experimental  design  with  mathematical 
methods  and  statistical  inferences  which  allows  the 
experimenter  to  empirically  analyze  the  system  of  interest. 

Since  the  relationship  between  the  response  and 
independent  variables  is  unknown,  the  first  step  is  to 
hypothesize  a  relationship  between  them.  In  many  cases, 
polynomial  models  are  used  as  the  approximating  function 
(10:399) . 

The  next  step  is  to  collect  the  data  based  on  an 
experimental  design.  The  method  of  least  squares  is  then 
applied  to  the  data  to  estimate  the  parameters  of  the 
functional  relationship.  This  regression  equation  can  then 
be  tested.  If  it  is  found  to  be  an  adequate  approximation  of 
the  true  functional  relationship,  the  experimenter  can  be 
confident  that  working  with  the  fitted  model  is 
representative  of  working  with  the  real  system. 

For  the  purposes  of  this  study,  it  was  hypothesized  that 
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the  functional  relationship  between  the  response  variables, 
aircraft  availability  and  sortie  generation  rate,  and  the 
independent  variables,  NTBF  and  MTTR,  was  a  second  order 
polynomial . 

The  experimental  design  chosen  was  a  second  order 
rotatable  central  composite  design.  This  design  requires 
five  levels  for  each  factor.  For  two  factors  the  design 

A 

consists  of  a  2  factorial  (coded  to  ±1  notation); 
augmented  by  4  axial  points  (^a,0),and  (0,jfea), 

where  a  =  (2*^)^/^  =  (2^)^/^  -  1.414; 
plus  n  center  points.  By  choosing  n  to  be  8,  the  central 
composite  design  sign  is  made  orthogonal  (10:462).  Thus, 
this  design  requires  4  runs  for  the  factorial,  plus  4  runs 
for  the  axial  points,  plus  8  runs  at  the  center  point  for  a 
total  of  16  runs.  Appendix  C  contains  a  layout  of  the 
design. 

Table  I  presents  the  five  levels  of  the  factors 
selected  for  the  design.  Looking  at  the  factorial  portion  of 
the  design,  a  change  of  50%  in  both  directions  results  in  a 
range  that  provides  a  wide  variation  for  determining  the 
effects  MTBF  and  MTTR  could  have  on  the  response  variables. 
However,  for  the  MTTR  data,  it  was  felt  that  the  low  level 
could  not  be  reduced  as  much  because  repair  work  requires 
some  minimum  amount  of  time.  Therefore,  only  a  10%  reduction 
of  the  Fairchild  Republic  data  was  selected  as  the  low  level 
of  MTTR.  This  still  results  in  a  range  that  provides 
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Table  I 

Levels  of  Factors 


MTBF 

MTTR 

Factorial  points 

High 

+  1 

150% 

150% 

Low 

-1 

50% 

90% 

Axial 

points 

High 

+1.414 

170.71% 

162.43% 

Low 

-1.414 

29.29% 

77.57% 

Center  points 

Center  0 

100% 

120% 

*  Asa 

percentage 

of  the  Fairchild 

Republic  data 

sufficient  variation  for  evaluating  the  effects  of  MTTR. 
Because  this  Is  a  central  composite  design,  the  center  point 
for  the  levels  of  MTTR  had  to  be  recomputed.  The  center 
point  turns  out  to  be  120%  of  the  Fairchild  Republic  data. 

The  sixteen  runs  are  then  performed  at  the  prescribed 
levels  of  the  factors.  After  collecting  the  data  from  the 
sixteen  runs,  regression  analysis  can  be  used  to  determine 
the  parameters  of  the  second  order  equation  representing  the 
data.  From  this  equation.  It  Is  possible  to  draw  response 


surfaces  which  describe  the  effects  of  the  Independent 
variables  on  the  response  variables. 


E  n  s  u  r  1  n 


ii 


Having  decided  which  factors  to  analyze  and  the  design 
which  would  provide  the  required  Information,  a  number  of 
questions  remained  to  be  answered  before  the  experimental 
design  could  be  performed.  There  were  five  particular 
questions  of  Interest. 

1.  How  long  of  a  warm-up  period  Is  required  to  avoid 
Initialization  bias? 

2.  How  much  time  Is  required  for  a  batch  mean 
observation? 

3.  Is  the  data  au tocor related? 

4.  How  many  batch  mean  observations  are  required? 

5.  What  Is  the  required  length  of  the  simulation? 


Warm-up  Period .  In  determining  how  long  of  a  warm-up 
period  was  required,  two  criteria  were  used.  The  first 
criteria  for  the  warm-up  period  was  to  determine  when,  over 
time.  Initialization  bias  no  longer  appeared  to  be  a  factor. 
An  Initial  run  was  made  for  7200  hours  (300  days).  The  plots 
of  aircraft  availability  and  sortie  generation  rate  versus 
time  were  observed  to  determine  when  the  values  reached  a 
s teady-s ta te .  By  2160  hours  the  values  of  both  aircraft 
availability  and  sortie  generation  rate  had  leveled  off 
Indicating  that  any  Initialization  bias  was  no  longer  a 
factor.  Thus,  the  suggested  warm-up  period  was  2160  hours. 

The  second  criteria  dealt  with  phase  Inspections.  Since 


as  part  of  the  simulation  Initialization  all  aircraft  were 


given  flight  hours  evenly  distributed  from  0  to  300  hours,  It 
was  felt  that  the  warm-up  period  should  allow  for  all 
aircraft  to  have  gone  through  the  phase  inspection  once.  For 
aircraft  flying  an  average  of  1.3  hours  per  sortie  and  an 
average  of  2.9  sorties  per  day,  it  would  require 
approximately  30  days  or  1920  hours  for  all  aircraft  to  go 
through  phase  inspection  once.  However,  since  the  first 
criteria  suggested  a  longer  warm-up  period,  it  was  concluded 
that  2160  hours  would  provide  a  valid  warm-up  period. 

Time  For  £  Ba  tch  Mean  Observation.  In  answering  the 
second  question  about  the  length  of  time  required  for  a  batch 
mean  observation,  the  batch  mean  observations  were  collected 
by  first  clearing  the  statistical  arrays  after  the  warm-up 
period  and  then  clearing  the  statistical  arrays  again  after 
every  set  amount  of  time  (e.g.  every  3  days).  The  batch  mean 
statistics  were  computed  by  taking  the  average  of  the 
statistics  over  the  time  that  they  were  collected. 

Determining  the  optimum  amount  of  time  required  for  these 
observations  involved  an  iterative  process  that  was  based  on 
a  number  of  factors:  whether  autocorrelation  was  present, 
the  calculated  number  of  batch  mean  observations  required, 
and  the  computer  time  which  would  be  needed  based  on  the 
computed  number  of  batch  mean  observations  required.  The 
discussion  on  autocorrelation  and  calculating  the  number  of 
observations  required  will  be  deferred  at  this  time  as  they 
are  covered  in  more  detail  in  questions  3  and  4  respectively. 


However,  separate  runs  were  made  which  cleared  the  statistics 
after  every  2,  3,  5,  10,  and  20  days.  Twenty  batch  mean 
observations  were  collected  for  the  2  and  3  day  periods,  and 
fifteen  batch  mean  observations  were  collected  for  the  5,  10, 
and  20  day  periods.  Based  on  tradeoffs  In  the  above  factors 
(autocorrelation,  number  of  observations,  and  computer  time). 
It  was  determined  that  the  batch  mean  statistics  collected 
over  10  days  (240  hours)  would  provide  the  accuracy  desired 
as  well  as  keeping  the  length  of  the  simulation  at  a 
manageable  level. 

Autocorrelation.  Detecting  autocorrelation  Is  Important 
since,  as  part  of  the  analysis  of  the  simulation,  a 
regression  model  is  used.  If  autocorrelation  is  present  and 
the  Ordinary  Least  Squares  (OLS)  estimators  are  used,  a 
number  of  problems  occur. 

1.  If  serial  correlation  is  allowed  for,  the 
estimators  will  be  Inefficient,  causing  the 
confidence  Intervals  to  be  unnecessarily  wide  and 
the  tests  of  significance  to  be  less  powerful. 

2.  If  autocorrelation  Is  Ignored  and  the  classical  OLS 
formulas  are  used,  the  usual  t  and  F  tests  of 
significance  are  no  longer  valid  and  if  used  would 
give  misleading  conclusions  about  the  significance 
of  the  estimated  regression  coefficients. 

3.  The  OLS  estimators  become  sensitive  to  sampling 
fluctuations  and  may  not  give  an  accurate  picture 
of  the  true  population  values  (9:226). 

Two  tests  were  used  to  detect  autocorrelation,  the 
Du r bln- Wa tson  d  test,  and  a  runs  test.  The  Durbin-Wa tson  d 
test  is  based  on  the  ratio  of  the  sum  of  the  squared 


differences  in  successive  residuals  to  the  Residual  Sum  of 
Squares  (9:235).  As  a  rule  of  thumb,  if  d,  the  test 
statistic,  is  found  to  be  2,  it  can  be  assumed  that  there  is 
no  first-order  autocorrelation,  either  positive  or  negative. 
The  closer  d  is  to  0,  the  greater  the  evidence  of  positive 
serial  correlation.  The  closer  d  is  to  4,  the  greater  the 
evidence  of  negative  serial  correlation  (9:237). 

The  BMDP  9R  program  was  used  to  calculate  the 
Durbin-Watson  statistic.  The  d  that  was  computed  was 
compared  against  critical  upper  and  lower  d  values  to 
determine  whether  autocorrelation  existed.  The  critical 
values  are  based  on  the  sample  size  and  number  of  explanatory 
variables . 

Using  the  BMDP  9R  program  on  observations  obtained  with 
a  batch  size  of  10  days,  the  computed  Durbin-Watson  statistic 
for  aircraft  availability  was  1.7349.  The  critical  values  at 
the  99%  level  of  confidence  with  n=15  and  k'=2  are 
dfower  "  0.70  and  =  1  .25  (9  :439  ).  Since  the  computed 

value  is  greater  than  d^ppg^  but  still  less  than  2,  we 
conclude  there  is  no  positive  serial  correlation.  The 
computed  Durbin-Watson  statistic  for  the  sortie  generation 
rate,  for  the  same  run,  was  2.4780.  The  critical  values  are 
the  same  as  above.  Since  the  computed  value  is  less  than 
2.75  (4  -  dyppgj,  =  4  -  1  .  25  =  2.75  )  but  is  greater  than  2,  we 
conclude  that  there  is  no  negative  serial  correlation. 

The  second  test  used  to  detect  autocorrelation  was  a 


runs  test.  Looking  at  the  residuals  for  the  same  15  batch 
mean  observations,  for  aircraft  availability  there  were  8 


runs  with  8  +'s  and  7  -'s.  The  critical  values  for  N1  =  8 
and  N2  =  7  Is  ^  4  or  2.  ^3  (9:440-441).  Since  8  Is  between 
these  values  we  conclude  that  there  Is  no  serial  correlation. 
Looking  at  the  residuals  for  sortie  generation  rate,  there 
are  10  runs  with  8  +'s  and  7  -'s.  The  critical  values  are 
the  same  as  above,  therefore,  we  again  conclude  that  there  Is 
no  serial  correlation  for  sortie  generation  rate.  Thus,  we 
are  relatively  confident  that  by  using  the  10-day  Interval  to 
determine  a  batch  mean,  autocorrelation  Is  not  present. 

Number  of  Ba  tch  Mean  Observa  tl ons  Required.  It  Is 
possible  to  obtain  a  value  of  an  output  variable  such  that  It 
estimates  the  true  population  value  within  some  accuracy 
criterion  with  a  high  degree  of  probability.  This  Is  done  by 
determining,  based  on  Initial  sample  values,  the  number  of 
observations  that  will  provide  the  desired  accuracy.  The 
number  of  batch  mean  observations  required  Is  determined  by 
the  following  formula  (6:427); 


where 


N  > 


N-l)(s) 


2 


(1) 


N  Is  the  number  of  observations  required, 

tg/2  f}_2  is  the  t-statlstlc  for  confidence  level  a/2 
’  and  N-1  degrees  of  freedom, 

s  Is  the  standard  deviation  of  the  sample,  and 

e  Is  the  half-width  of  the  confidence  Interval. 


A  confidence  level  of  95%  (a=0.05)  and  half  width  of 


*0.005  for  both  aircraft  availability  and  sortie  generation 
rate  were  used.  The  half  width  for  aircraft  availability 
yields  an  estimate  to  within  1%.  The  half  width  for  sortie 
generation  rate  results  In  a  rate  to  the  nearest  0.01.  These 
half  widths  were  considered  reasonable  for  this  study  and 
were  used  to  determine  the  number  of  observations  required. 

Table  II  contains  the  calculations  for  aircraft 
availability  and  sortie  generation  rate  using  the  10-day 
batch  size  and  13  batch  means. 

Table  II 


Calculations  for  Number  of  Observations 


Aircraft  Availability 

Sortie  Generation  Rate 

n=15,  a  =0.05 

n=15,  a  =0.05 

*^0.025,14  “  2.145 

*^0.025,14  “  2.145 

X  =  0.91302 

X  =  2.92833 

S  =  0.00329 

S  =  0.00740 

e  =  0.005 

e  =  0.005 

N  =r(  2.145)(0.00329r|2 

N  =  R  2. 145)  (0.00740)“]  2 

L  (0.005)  J 

L  (0.005)  J 

N  =  2* 

N  *  11* 

*  Rounded  up  to  the 

nearest  Integer 

sample,  only  11  batch  mean  observations  are  needed.  This  Is 
the  governing  factor  for  the  number  of  observatlor = .  Using 
the  10-day  period  with  11  observations  requires  that  the 
simulation  be  run  for  2640  (240  x  11)  hours. 

Simulation  Length.  The  length  of  the  simulation  is 
determined  by  the  amount  of  warm-up  time  required  plus  the 
amount  of  time  required  for  data  collection.  For  this  study 
the  warm-up  period  was  determined  to  be  2160  hours  and  the 
amount  of  time  required  for  data  collection  was  determined  to 
be  2640  hours.  Therefore,  the  total  length  of  the  simulation 
should  be  4800  hours  (2160  +  2640)  or  200  days. 


IV.  Results 


The  last  chapter  presented  the  design  chosen  to 
Investigate  the  effects  of  MTTR  and  MTBF  on  aircraft 
availability  and  sortie  generation  rate.  In  addition,  a 
number  of  questions  relevant  to  running  the  simulation  were 
discussed.  This  chapter  discusses  the  results  of  the 
analysis  performed  on  the  model. 

As  was  discussed,  a  rotatable  central  composite  design 
was  chosen  to  obtain  data  points  for  the  regression  analysis. 
This  design  required  16  runs  to  be  performed.  The  coded 
levels  of  the  factors  as  well  as  the  results  for  each  run  are 
detailed  In  Appendix  C. 

The  BMDP  9R  stepwise  regression  program  was  run  on  both 
the  aircraft  availability  and  sortie  generation  rate  data 
using  a  second  order  equation  of  the  form: 

Y  =  Bq  +  BjXj  +  B2X2  +  B12X1X2  +  +  B22X2^  (2) 

where  X^  and  X2  are  MTTR  and  MTBF  respectively. 

The  'Best'  subset  of  the  five  variables  was  determined 
using  the  Mallows'  Cp  criterion.  The  Mallows'  Cp  criterion 
attempts  to  Identify  the  subset  of  variables  that  has  the 
smallest  total  mean  squared  error.  When  the  Cp  value  for 
this  subset  Is  also  near  p  (the  number  of  parameters  In  the 
model),  the  bias  of  the  regression  model  Is  small  (11:426-427) 

Alrcraf  t  Avallablll tv  Results 

When  the  ANOVA  (Analysis  of  Variance)  for  aircraft 


47 


i 

f 

r 

1 

*  7 

;  availability  based  on  a  2  factorial  was  performed,  both 

factors,  MTTR  and  MTBF,  as  well  as  their  Interaction  term 
were  significant.  The  ANOVA  results  are  presented  In 
Table  III. 


Table  III 


ANOVA  Table  for  Aircraft  Availability 


Source 

SS 

df 

MS 

F 

MTTR 

0.0188701 

1 

0.0188701 

409.01* 

MTBF 

0.0705441 

1 

0.0705441 

1529.03* 

Interac  tlon 

0.0068650 

1 

0.0068650 

148.80 

Error 

0.0018455 

40 

0.0000461 

0.1030216 

43 

*  Significant 

a  t 

a  =  1%  level 

However,  In  running  BMDP  9R  on  the  16  data  points 
obtained  from  the  central  composite  design,  the  'Best'  subset 
was  only  In  terras  of  the  MTBF  variable  (the  Intercept,  MTBF, 
and  MTBF  squared).  For  this  subset.  Mallows'  Cp  was  2.92  for 
the  three  parameters  In  the  model  (p  =  3).  This  would 
Indicate  that  there  Is  little  bias  in  the  model.  The 
regression  Indicated  that  aircraft  availability  can  be 
explained  by  the  following  equation  (See  Table  D.l  of 
Appendix  D) : 

AA  =-  0.917869  +  0 . 07  3066  2  (  X2  )  -  0 . 0  547  96  (  Xj  ^  )  (3) 


whe  re 


AA  Is  aircraft  availability,  and 
Xo  Is  the  coded  level  of  MTBF. 


This  equation  explains  approximately  81%  of  the  variation  In 
aircraft  availability.  Figure  4  provides  a  graphical 
representation  of  this  equation.  As  might  be  expected,  when 
MTBF  Increases  (l.e.  when  there  Is  more  time  between 
failures)  the  availability  Increases. 


AIRCRAFT 

AVAILABILITY 


LEVELS  OF  MTBF 


*  As  a  percentage  of  the  Fairchild  Republic  data 


Figure  4.  Aircraft  Availability  for  Various  Levels  of  MTBF 


Of  particular  Interest  on  this  graph  Is  the  level  of 
MTBF  which  results  In  83%  aircraft  availability.  This  level 
Is  the  minimum  contract  requirement  for  availability  (4:5). 


The  levels  of  MTBF  that  result  In  less  than  83%  availability 


can  be  determined  from  Figure  4.  An  availability  of  83% 
corresponds  to  a  decrease  in  the  MTBF  level  of  approximately 
38%. 

Sortie  Generation  Rate  Results 

The  ANOVA  for  SGR  based  on  the  2^  factorial,  Indicated 
that  both  factors,  MTTR  and  MTBF,  were  significant  at  the  95% 
level  but  that  the  Interaction  term  was  not  significant.  The 
ANOVA  results  are  contained  Table  IV. 


Table  IV 

ANOVA  Table  for  Sortie  Generation  Rate 


Source 

SS 

df 

MS 

F 

MTTR 

0.000596459 

1 

0.000596459 

8.57* 

MTBF 

0.000349455 

1 

0.000349455 

5.02 

Interac  tlon 

0.000160366 

1 

0.000160366 

2.31 

Error 

0.002782899 

40 

0.000069572 

0.003889179 

43 

*  Significant 

a  t 

a  =  5%  level 

However,  after  running  BMDP  9R  on  the  data,  the  'Best' 
subset  was  again  only  in  terms  of  the  MTBF  variable  (the 
intercept,  MTBF,  and  MTBF  squared).  For  this  subset. 

Mallows'  Cp  was  1.04  for  the  three  parameters  indicating  that 
there  was  slightly  more  bias  in  this  model  than  there  was  for 
aircraft  availability.  The  regression  Indicated  that  sortie 
generation  rate  can  be  explained  by  the  following  equation 


Of  particular  Interest  on  this  graph  Is  the  fact  the 
predicted  sortie  generation  rate  never  decreases  to  2.2.  The 
T-46A  Mill  be  flown  at  an  average  rate  of  60  hours  per  month 
per  aircraft  (4:6).  This  utilization  rate  Is  equivalent  to  a 
2.2  sortie  generation  rate. 

The  main  reason  the  regression  equation  (4)  explains 
only  51%  of  the  variation  Is  due  to  the  manner  In  which 
flight  line  operations  are  modeled.  The  simulation  Is  set  up 
to  allow  all  82  aircraft  to  fly  with  the  major  restriction  on 
the  sortie  rate  being  the  stratified  takeoff  (l.e.  allowing 
only  one  aircraft  to  takeoff  every  three  minutes).  Due  to 
this  limitation  there  Is  almost  always  an  aircraft  waiting  to 
takeoff,  even  when  the  rate  at  which  aircraft  fall  Is 
Increased.  Because  the  system  Is  not  stressed,  the  sortie 
generation  rate  tends  toward  the  maximum  allowable  rate  per 
day.  If  only  a  portion  of  the  aircraft  were  allowed  to  fly 
all  the  sorties  for  a  given  day,  the  MTBF  would  possibly 
account  for  a  higher  percent  of  the  variation  In  the 
resultant  sortie  generation  rate  equation. 

As  was  mentioned,  the  regression  equation  for  both 
aircraft  availability  and  sortie  generation  rate  did  not 
Include  a  MTTR  factor.  This  is  not  to  Imply  that  MTTR  is  not 
a  significant  factor.  The  ANOVA  Indicated  that  It  was. 
However,  In  performing  the  regression  analysis,  once  MTBF 
entered  the  equation,  MTTR  was  unable  to  explain  any  of  the 
remaining  variance.  Even  when  MTTR  is  forced  Into  the 


equation,  the  coefficients  for  MTTR  are  statistically  not 
different  from  zero.  Therefore,  MTTR  did  not  show  up  in 
either  best  fit  equation  because  it  was  unable  to 
substantially  contribute  to  the  explanation  of  the  remaining 
variation  once  MTBF  entered  the  equation. 


V .  Concl u s 1 ons  and  Recommendations 

Conclusions 

Based  on  the  verification  and  validation  efforts  for 
this  model,  our  conclusion  Is  that  this  Is  a  valid  model  for 
predicting  the  availability  and  sortie  generation  rate  of  the 
T-46A.  Due  to  the  manner  In  which  this  model  has  been 
constructed,  It  Is  a  flexible  model  that  can  be  easily 
adapted  to  different  aircraft. 

From  the  analysis  outlined  In  Chapter  IV  the  research 
questions  posed  In  Chapter  I  can  now  be  answered. 

1)  What  Is  the  availability  of  the  T-46A  In  the  UPT 
system?  The  average  aircraft  availability  predicted  by  the 
model  is  0.921.  This  estimate  uses  Fairchild  Republic 
predicted  failure  rates  as  opposed  to  allocated  rates.  A  95% 
confidence  Interval  for  this  prediction  is  [0.917,  0.923). 

2)  What  sortie  generation  rate  can  be  achieved  by  the 
T-46A?  The  average  sortie  generation  rate  predicted  by  the 
model  Is  2.926.  This  again  Is  using  the  Fairchild  Republic 
predicted  failure  rates.  A  95%  confidence  Interval  for  the 
sortie  generation  rate  prediction  is  [2.923,  2.932). 

The  estimated  value  of  the  sortie  generation  rate  Is  the 
maximum  that  can  occur  under  the  stratified  takeoff 
restrictions.  Only  In  the  most  pessimistic  case  does  the 
sortie  generation  rate  show  a  modest  drop.  When  the  failure 
rate  Is  Increased  by  more  than  70%,  the  sortie  generation 


rate  drops  to  2.715.  However,  this  rate  Is  still  much  above 
the  projected  utilization  rate  of  60  hours  per  month  per 
alrcraf  t. 

In  contrast  to  the  sortie  generation  rate,  the  estimates 
of  aircraft  availability  did  not  meet  contract  requirements 
at  all  levels  of  the  experiment.  Analysis  shows  that  If  the 
failure  rates  Increase  by  38%  over  the  contractors  estimate, 
the  T-46A  will  be  unable  to  achieve  an  average  availability 
of  83%. 

In  summary,  the  factor  that  has  the  most  Influence  on 
both  availability  and  sortie  generation  rate  Is  the  mean  time 
between  failures.  Although  analysis  of  variance  shows  mean 
time  to  repair  to  be  a  significant  factor,  It  does  not  enter 
Into  either  regression  equation  used  to  predict  availability 
or  sortie  generation  rate. 

Recommendations 

There  are  three  areas  that  suggest  further  effort.  The 
first  area  Is  scheduled  maintenance  and  how  It  Is  accounted 
for  In  the  model.  The  second  area  Is  manpower  and  how  to 
adjust  levels  to  reflect  only  T-46A  work  and  not  T-38  work, 
training,  and  other  utilization  factors.  The  final  area  Is 
the  Inclusion  of  spares  and  off-alrcraft  maintenance  In  the 
mode  1 . 

Scheduled  maintenance  Is  accounted  for  In  this  model  by 
delaying  the  aircraft  once  every  300  flying  hours  for  phase 
Inspection.  All  maintenance  manhours  associated  with 


scheduled  maintenance  actions.  Including  those  not  at  300 
hour  Intervals,  are  accumulated  when  the  aircraft  enters  the 
phase  Inspection  network.  Accounting  for  these  actions  as 
they  are  scheduled  to  occur  would  give  a  more  accurate 
representation  of  the  UPT  system.  A  flying  hour  counter  for 
each  scheduled  maintenance  action  would  be  necessary.  This 
would  require  fifteen  counters  In  addition  to  the  one 
currently  used  to  send  an  aircraft  to  phase. 

Another  recommendation  Involves  manpower.  We  were 
unable  to  accurately  breakdown  the  manpower  levels  given  In 
the  ATC  LCOM  Final  Report.  The  levels  contained  there 
reflect  not  only  the  on-alrcraft  maintenance  but  also 
off-alrcraft  and  T-38  maintenance,  and  factors  for  training, 
leave,  sick  days,  etc.  The  manpower  normally  used  for  these 
other  activities  Is  used  for  on-aircraft  maintenance  In  our 
model  of  the  T-46A.  Thus,  the  manpower  levels  were  not 
constraining  factors  as  expected.  That  manpower  is  not  at  a 
constraining  level  may  be  seen  by  comparing  the  estimate  of 
availability  given  above  (0.921)  to  an  estimate  obtained  with 
an  unlimited  manpower  level.  Using  a  level  of  10,000 
personnel  for  each  specialty  code,  the  estimate  for 
availability  Is  0.922  with  a  95%  confidence  interval  of 
[0.920,  0.924].  This  estimate  Is  not  significantly  different 
from  the  previous  estimate.  Further  research  into 
determining  appropriate  manpower  levels  Is  recommended. 

The  final  recommendation  deals  with  spares  and 


A  .■  A  A  '  A  'i^V^  W^i.^  -u*r.- '■^;  i.-  rf-  ir-7«-.i  •  -  *-. 


off-aircraft  maintenance.  This  effort  at  modeling  the  T-46A 
was  directed  at  only  on-aircraft  maintenance.  In  other 
words,  the  model  assumes  that  spares  are  always  available  and 
that  off-aircraft  maintenance  is  independent  from  on-aircraft 
maintenance.  The  effects  of  spares  levels  and  off-aircraft 
maintenance  may  be  of  interest.  The  model  could  be  expanded 
to  include  these  considerations. 
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T-46A  Model 


This  appendix  contains  the  SLAM  and  FORTRAN  codes  of  the 
T-46A  model  developed  In  this  study.  The  first  section  lists 
the  SLAM  code.  The  second  section  is  an  explanation  of  the 
allocation  subroutine  contained  in  the  FORTRAN  code.  The 
final  section  lists  the  FORTRAN  code. 

SLAM  Code 


GEN, FOLEY  HAGER, T46  MODEL ,10/10/85 ,1  ,N,N; 

;  ATTRIBUTES  USED 

;  1—  STORES  TNOW  TO  COLLECT  DOWN  TIME 

;  2—  TAIL  NUMBER 

;  3—  NUMBER  OF  FLYING  HOURS 

;  4—  STORES  NUMBER  OF  FAILURES 

;  5—  STORES  REPAIR  TIME 

;  6“  STORES  REPAIR  CREW  CODE 

;  7—  STORES  NUMBER  OF  REPAIR  PERSONS  USED 

;  8--  STORES  LENGTH  OF  PREFLIGHT  AND  THRUFLIGHT 

;  9—  STORES  WUC  OF  FAILURE  WHEN  THERE  IS  A  RMA 

;  10—  STORES  WHICH  RMA  IS  BEING  WORKED  ON 

;  11—  STORES  WHICH  FAILURE  NUMBER  THE  FAILURE  IS 

;  12—  STORES  TIME  OF  LAST  PREFLIGHT 


FILES/QUEUES  USED 

1—  AIRCRAFT  WAIT  FOR  F43!.  TO  PERFORM  PREFLIGHT  INSPECTION 

2—  AIRCRAFT  WAIT  FOR  DAYLIGHT 

3— 

4— 

5—  AIRCRAFT  WITH  FAILURES  WAIT  FOR  REPAIR  CREWS 

6—  AIRCRAFT  WAIT  FOR  PHASE  DOCKS 

7—  AIRCRAFT  WAIT  FOR  STRATIFIED  TAKE  OFFS 

8—  AIRCRAFT  PLACED  IN  THIS  QUEUE  TO  BEGIN  FLIGHT  SEQUENCE 

9—  AIRCRAFT  WAIT  FOR  F431  TO  PERFORM  POSTFLIGHT  SERVICING 

10—  AIRCRAFT  WAIT  FOR  SCHEDULING  TO  OCCUR 

11— 

12—  TAKE  OFF  CREATION  ENTITY  WAITS  HERE  FOR  DAYLIGHT 

13—  SLAM  CALENDER  OF  EVENTS 
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LIMITS, 12, 12, 200; 

INTLC,XX(1)=0; 

INTLC,XX(2)=0; 

INTLC,XX(3)=0; 

INTLC,XX(4)=0; 

INTLC,XX(5)=0; 

INTLC,XX(6)=0; 

INTLC,XX(7)=0; 

INTLC,XX(8)=0; 

INTLC,XX(9)=0; 

INTLC,XX(10)*0; 

INTLC,XX(11)=0; 

INTLC,XX(12)=0; 

INTLC,XX(13)=0; 

INTLC,XX(14)=0; 

INTLC,XX(99)=0; 


NUMBER  OF  AIRCRAFT  FULLY  MISSION  CAPABLE 

TOTAL  MAINTENANCE  MANHOURS 

TOTAL  NUMBER  OF  FLIGHTS 

SORTIE  GENERATION  RATE 

COUNTER  TO  INITIALIZE  FLT  HOURS 

SORTIE  DURATION 

MMH/S 

TOTAL  DOWNTIME 
DOWN  TIME  PER  SORTIE 

ACCUMULATE  REPAIR  HOURS  TO  COMPUTE  MTTR 
COUNTS  NUMBER  OF  SWITCHES  IN  ALLOC  SUBROUTINE 
AVERAGE  DOWN  TIME 

AVG  NUMBER  OF  AIRCRAFT  FULLY  MISSION  CAPABLE 
NUMBER  OF  FLIGHTS  SINCE  LAST  STATS  COLLECTION 
TOTAL  NUMBER  OF  AIRCRAFT 


TIMST,XX(1),FMC  AIRCRAFT; 

TIMST,XX(4), SORTIE  RATE; 

TIMST,XX03)  .AVERAGE  AC  AVL;  XX(13)  IS  SET  =  TTAVG(l)  [XX(1)]  IN  EVENT  3 


RECORD, TNOW, TIME, 0,T, 240, 240;  SAMPLING  EVERY  240  HOURS  IF  THIS  CHANGES 
VAR, XX( 13), A, AVG  AVAIL  ACFT;  MUST  CHANGE  DIVISOR  IN  NODE  SGR 
VAR,XX(4) ,S,SORTIE  RATE; 

PRIORITY/ 5, LIFO; 

;  TIME  UNIT  IS  ONE  HOUR 

NETWORK; 


IDENTIFY  RESOURCE  CONSTRAINTS 


CREW  RESOURCES 


RESOURCE/R431(0),5; 
RESOURCE/F431(0) ,1 ,5 ,9; 
RESOURCE/P431(0),5; 

RES0URCE/S4270(0),5; 
RES0URCE/S4274(0),5; 
RESOURCE/S4275(0) ,5; 
RESOURCE/S4271(0),5; 


T-37  REPAIR  AND  RECLAMATION 

T-37  FLM 

T-37  INSPECTION 

MACHINE 

METALS  PROCESSING 
STRUCTURAL  REPAIR 
CORROSION  CONTROL 


RESOURCE/S426(0) ,5 
RESOURCE/ A426(0), 5 
RESOURCE/T426(0)  ,5 
RESOURCE/F426(0),5 


T-37  JET  ENGINE  SHOP 

T-37/38  ACCESSORY  REPAIR 
T-37/38  TEST  CELL 
T-37  FLIGHT  LINE  SUPPORT  UNIT 


RESOURCE/ W431(0),5; 


T-37/38  WHEEL  AND  TIRE 


RESOURCE/S4233(0) ,5;  T-37/38  FUEL  SYSTEMS 

RESOURCE/S4230(0),5;  T-37/38  ELECTRICAL  SYSTEMS 

RESOURCE/S4234(0) ,5;  T-37/38  PNEUDRAULICS  SYSTEMS 

RESOURCE/S4231(0),5;  T-37/38  ENVIRONMENTAL  SYSTEMS 

RESOURCE/S4232(0) ,5;  T-37/38  EGRESS  SYSTEMS 

RESOURCE/S3280(0),5;  T-37/38  RADIO  AND  RADAR  REPAIR 

RESOURCE/S3281(0) ,5;  T-37/33  RADIO  AND  RADAR  REPAIR 

RESOURCE/S325(0),5;  T-37/38  AUTO  FLIGHT  CONTROL 

;  OTHER  RESOURCES 

RESOURCE /DOCK (4), 6; 

GATE/ DAY , CLOSED , 2 , 1 2 ; 

GATE/ STRAT , CLOSED , 7 ; 

GATE/SCHEDULE , CLOSED , 10 ; 

;  **********  main  network  ********** 

CREATE, 0,, 1,82;  GENERATE  82  AIRCRAFT 

ASSIGN, XX(1)=XX(1)+1; 

ASSIGN,ATRIB(2)=XX(1);  ASSIGN  TAIL  NUMBERS 

ASSIGN, XX(99)=XX(99)+1;  COUNT  AIRCRAFT 

ASSIGN, ATRIB(3)=XX(5); 

ASSIGN, XX(5)=XX(5)+3. 7;  INITIALIZE  FLT  HOURS 

j *********************************************************************** 


» 

;*  MODEL  SEGMENT  I  ************  SORTIE  GENERATION  *****************  * 


SCH  AWAIT( 10), SCHEDULE; 

APRE  AWAIT(1),F431; 

SPFL  ASSIGN, ATRIB(8)  =  RLOGN( 1 . 414219 , . 
A^^rll  ASSIGN, XX(2)  =  XX(2)  +  ATRIB(8); 

ACT,ATRIB(8); 

FPRF  FREE,F431; 

ASSIGN, ATRIB( 12)  =  TNOW; 

AWDAY  A WAIT (2), DAY; 

SRTTO  AWAIT(7),STRAT; 

TFLY  QUEUE(8); 

ACT(82)/1 ,.2; 

SLEN  ASSIGN,XX(6)=RNORM(l .3 , .13 ,5) , 
ATRIB(3)=ATRIB(3)+XX(6); 
ACT/2,XX(6); 

SSTS  ASSIGN, ATRIB(l)=TNOW, 

XX(14)  =  XX(14)  +  1, 
XX(3)=XX(3)+1; 

ACT, .2; 


WAIT  FOR  CREW  CHIEF 
551545,4);  SET  PREFLIGHT  LENGTH 
ACCUMULATE  MMH 
PERFORM  PREFLIGHT 
RELEASE  CREW  CHIEF 
STORE  PREFLIGHT  TIME 
WAIT  FOR  DAYLIGHT 


TAXI  FOR  LAUNCH 


SET  SORTIE  LENGTH 
FLY  MISSION 


COUNT  SORTIES 
AFTER  FLIGHT  TAXI 


PSTFL  AWAIT(9),F431; 

STFL  ASSIGN, ATRIB(8)  =  RLOGN( .435 16965,6) 
AM2  ASSIGN, XX(2)  =  XX(2)  +  1 .326*ATRIB{8) ; 

ACT,ATRIB(8); 

FPSTF  FREE,F431; 

MMHPS  ASSIGN, XX(7)=XX(2)/XX(3); 

DNTM  ASSIGN, XX(9)  =  XX(8)/XX(3); 

SGR  ASSIGN, XX(4)=XX(14)/XX(99)/ 10; 

F?PH2  G00N,1; 

ACT,,ATRIB(3).GE.300.0,PHAZ; 

ACT,,, FAIL; 

CONT  G00N,1; 

ACT , , NNGAT( DAY ) . EQ , 0 , SRTT ; 

ACT, , ,SCH; 


AWAIT  CREW  CHIEF 
;  SET  THRUFLIGHT  LENGTH 

ACCUMULATE  MMH 
PERFORM  THRUFLIGHT 

SET  MMH/SORTIE 
SET  DOWNTIME  PER  SORTIE 
SET  SORTIE  GENERATION  RATE 

IS  PHASE  INSPECTION  DUE? 

IF  NOT,  CHECK  FOR  FAILURE 

FLY  AGAIN  IF  STILL  DAY 
OR,  SCHEDULE  FOR  NEXT  DAY 


jHr*ir'jHr********ir**************-jHtr************'<HM>;*******************i>r***^Hr 


MODEL  SEGMENT  II 


****  CHECK  FOR  FAILURES  *** 


FAIL 

F?REP 


DETC 

FISTS 

REP 

MTTR 

AMH3 

RELCR 

MLTR? 


F2STS 


EVENT , 1 ; 

G00N,1; 

ACT,, ATRIB(4).EQ. 0.0, CONT; 

ACT; 

EVENT, 4;  DETERMINE  WHICH  CREW 

ASSIGN, XX(l)=>XX(l)-l,ATRIB<l)=TNOW; 
AUAIT(5),ALLOC(2); 

ASSIGN, XX(10)=XX(10)+ATRIB(5); 

ASSIGN, XX(2)=XX(2)+ATRIB(5)*ATRIB(7); 
ACT,ATRIB(5); 

EVENT, 2; 

GOON , 1 ; 

ACT,,ATRIB(4).GT.0.0,REP; 

ACT; 

ASSIGN,XX(1)=XX(1)+1 , 

XX( 12 )=TNOW-ATRIB( 1 ) , 
XX(8)=XX{8)+XX(12); 

GOON , 1 ; 

ACT , ,TNOW-ATRIB( 12 ) . LT. 24. 0 ,CONT; 
ACT,,,SCH; 


BACK  TO  FLY  IF  NO  REPAIRS 

IS  NEEDED  TO  REPAIR  FAILURE 

WAIT  FOR  REPAIR  CREWS 
COLLECT  STATS  ON  MTTR 
ACCUMULATE  MMH 
REPAIR  ACTIVITY 
RELEASE  REPAIR  CREWS 


COLLECT  DOWNTIME  STATS 
BACK  TO  CONT  IF  PREFLIGHT 
IS  NOT  NEEDED 
TO  SCHEDULE  IF  ONE  IS 


MODEL  SEGMENT  III 


******  phase  INSPECTION  ****** 


*********************************************************************** 


PHAZ  AWAIT(6),DOCK; 

PISTS  ASSIGN, XX(1)=XX(1)-1; 

AMH4  ASSIGN ,XX( 2)=XX( 2)+97 . 79845 ; 
RSa  ASSIGN,  ATRIB(  3  )=0; 


ACCUMULATE  ALL  SCHEDULED  MMH 
RESET  PHASE  CLOCK 


%V. 


. 
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ACT/2,72; 

FDOCK  FREE, DOCK; 

P2STS  ASSIGN, XX(1)=XX(1)+1; 
ACT, , ,SCH; 


3-DAY  PHASE  DURATION 


TO  SCHEDULE 


*  * 

*  MODEL  SEGMENT  IV  ******  DAY/NIGHT  ******  * 

*  * 


CREATE; 

ACT,1; 

EVENT, 3; 

ACT, 7; 

SRISE  OPEN, DAY; 

ACT, 12; 

CSCH  CLOSE, SCHEDULE; 
NIGHT  CLOSE, DAY; 

“ACT,  4; 

SCHPL  EVENT, 3; 

ACT, 8,, SRISE; 


12  HRS  OF  DAYLIGHT 

CLOSE  GATE  TO  COLLECT  AIRCRAFT 

NIGHT 

SCHEDULE  WHICH  PLANES  FLY  THE  NEXT  DAY 


MODEL  SEGMENT  V 


******  WORK  SHIFTS  ****** 


* 
* 

;*  * 
•  *********************************************************************** 

CREATE,, 8; 

DSHF  EVENT, 6; 

ACT, 8; 

SSHF  EVENT, 7; 

ACT, 8; 

NSHF  EVENT, 8; 

ACT, 8,, DSHF; 

*********************************************************************** 
*  * 

*  MODEL  SEGMENT  VI  ******  CREATE  STRATIFIED  TAKE  OFFS  ******  * 

*  * 


DAY  SHIFT  COMES  ON  DUTY 
SWING  SHIFT  COMES  ON  DUTY 
NIGHT  SHIFT  COMES  ON  DUTY 


CREATE; 

TOAD  AWAIT(12),DAY; 
TO  EVENT, 5; 

ACT,. 05,, TOAD; 


ALLOW  ONE  PLANE  TO  TAKE  OFF 
DELAY  NEXT  TAKE  OFF  FOR  3  MINUTES 


MODEL  SEGMENT  VII 


******  CLEAR  STATS  ****** 


* 

* 

* 


CREATE,, 6.0; 

CLER  ASSIGN, XX( 14)  =  0; 
ACT,240,,CLER; 


ENDNETWORK; 
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INIT, 0,4800; 

MONTR, CLEAR, 0.001 ,240.0; 

FIN; 

Explanation  of  Alloc  Subr ou  tine 

The  allocation  subroutine  Is  used  to  allocate  crew 
resources  for  the  repair  of  aircraft  failures.  In  order  to 
understand  the  structure  of  this  subroutine  a  description  of 
how  SLAM  calls  this  subroutine  must  be  given  first.  SLAM 
calls  the  allocation  subroutine  when  either  of  two  events 
happen.  The  first  event  Is  the  arrival  of  an  entity  at  an 
await  node  tied  to  the  allocation  subroutine.  The  second 
event  Is  a  change  in  the  quantities  of  any  of  the  resources 
used  by  the  allocation  subroutine.  It  Is  also  Important  to 
note  that  SLAM  only  checks  one  entity  when  the  allocation 
subroutine  Is  called.  Which  entity  is  checked  depends  on  the 
queue  discipline  for  that  file. 

In  the  model,  aircraft  can  arrive  at  the  node  allocating 
crew  resources  from  two  locations.  The  first  location  Is  the 
SLAM  node  labeled  FISTS.  In  this  case  It  has  just  been 
determined  that  the  aircraft  has  a  failure.  The  second 
location  Is  the  SLAM  node  MLTR?.  Here  the  aircraft  has  more 
maintenance  actions  to  be  accomplished.  This  could  be  the 
result  of  the  aircraft  having  multiple  failures  or  the 
maintenance  action  just  completed  Is  a  required  maintenance 


action. 


Thus,  aircraft  can  arrive  at  the  allocation  node  in  one 
of  two  states.  There  is  either  no  work  in  progress  (a  state 
called  NOWIP)  or  the  aircraft  has  just  finished  a  required 
maintenance  action  (state  RMAWIP).  When  the  allocation 
subroutine  is  called  a  check  is  made  to  determine  which  state 
the  aircraft  is  in.  Based  on  the  state,  a  call  is  made  to 
either  the  NOWIP  or  RMAWIP  subroutine.  Figure  A.l  shows  a 
flow  chart  of  the  allocation  subroutine.  Both  the  NOWIP  and 
RMAWIP  subroutines  provide  the  allocation  subroutine  with  the 
crew  code  and  size  and  the  repair  time  for  the  failure 
selected  to  be  fixed.  If  there  are  no  crews  available  to  fix 
a  failure,  the  crew  size  is  zero  and  the  crew  code  and  repair 
time  are  meaningless.  Next,  the  allocation  subroutine  checks 
to  see  if  a  crew  is  available.  If  a  crew  is  available,  the 
subroutine  seizes  the  appropriate  number  of  the  crew 
available  to  fix  the  failure.  When  there  are  no  crews 
available,  the  allocation  subroutine  checks  to  see  if  there 
are  more  aircraft  waiting  for  repairs  that  need  to  be 
checked.  If  there  are,  the  order  of  the  aircraft  in  the  file 
are  switched  to  allow  another  aircraft  to  be  checked. 

Subroutines  NOWIP,  RMAWIP  and  the  subroutines  called 
by  them  (except  REDUCE)  are  presented  only  as  flowcharts  in 
Figures  A. 2  through  A. 5. 


Are  required  maintenance 
actions  being  done 
to  this  aircraft? 


Yes 


No 


Call  NOWIP 


Call  RMAWIP 


Is  a  crew  available? 


Yes 


1 


No 


Have  all  aircraft 
waiting  for  crews 
been  checked? 


Yes 


No 

f 

Switch  order 
of  alrcraf  t 


Seize  resources 

I 

store  repair  time, 
crew  code,  and 
crew  size  In 
attributes  5-7 


Re  turn 


Does  this  failure  have 
required  maintenance 
actions? 


Yes 


Is  next 
activity  code 
1? 


Get  crew  data  for 
required 

maintenance  action 


Get  crew  data  for 
possible 

maintenance  action 


Are  there  enough 
repair  personnel 
available  for 
repair? 


Set  ICSZND  =  0 


Is  next 
activity  code 
1? 


Set 

a  t  tr i bu  te  s 
9-11  to  0 


Store  required 
main  tenance 
data  in 
attributes 
9-11 


Set  NEED.ISIZE 
and  RTIME 


Return 


Figure  A. 4  Subroutine 

AVAIL(NPLANE  ,NACT, IFLNUM.MA  ,NEXTRA  ,NEED ,ISIZE,NEED) 
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Is  there  more  than 
one  failure  to  check? 


Crew  code  =  NEED 


Crew  Code  = 
Code  of  failure  I 


Is  crew 

going 

to  be  i 

ised? 

Has  the  aircraft 
been  waiting 
for  a  crew? 


Begin  collection 
of  waiting 
statistics  for 
this  crew 


Has  the  Yes 

aircraft - 1 

been  waiting  ▼ 

for  the  crew?  End  collection 

of  waiting 
statistics  for 
No  this  crew 


Are  there 
more  failures? 


Count  number 
of  uses 


Re  turn 


Figure  A. 5  Subroutine  CREWSTATS( NPLANE ,NCHECK ,NEED  ,  ISIZE  ) 
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FORTRAN  Code 


C^r^rk^t^rk^e^c^rk^e^rk^[^e^rk^t^e^r^e^rk^eirk^rk^t^c*^^^e■k^He^c^c/rk^e^rk^e^t^rk*irk**ie^^^e^rk^rk*1e1e^e^eie**^^^rk 

* 

MAIN  PROGRAM  * 

* 

•*★*********★*****★*★*****★★****★******* 
C  THIS  PROGRAM  DEFINES  ALL  GLOBAL  VARIABLES 

C  THE  NAME  OF  THE  USER  VARIABLES  ARE  EXPLAINED  BENEATH  EACH  COMMON 
C  DECLARATION.  ALL  NAMES  ARE  IMPLICITLY  TYPED. 

DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

COMMON/SCOMl/  ATRIB(IOO)  ,DD(  100)  ,DDL(  100)  ,DTNOW, II ,MFA .MSTOP  ,NaNR 
1 ,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOMl/  ITOTFL .NCREWS ,NW0C 
C  ITOTFL  —  COUNTS  TOTAL  NUMBER  OF  FAILURES 
C  NCREWS  --  NUMBER  OF  DIFFERENT  SPECIALISTS  MODELED 
C  NWUC  --  NUMBER  OF  DIFFERENT  WUCS  MODELED 
COMMON/UCOM2/  ISHFTS(20,3) 

C  ISHFTS  —  CONTAINS  THE  CHANGE  IN  MANPOWER  FROM  ONE  SHIFT  TO  THE  NEXT 
COMMON/UCOM3/  NFMAT( 74 ) ,NPMA( 74) ,NRMA( 74) 

C  NFMAT  --  COUNTS  FAILURES  BY  WUC 

C  NPMA  --  NUMBER  OF  POSSIBLE  MAINTENANCE  ACTIONS  BY  WUC 
C  NRMA  —  NUMBER  OF  REQUIRED  MAINTENANCE  ACTIONS  BY  WUC 
C0MM0N/UC0M4/  NCCRA( 74,2) ,NCSRA( 74,2) 

C  NCCRA  --  CREW  CODES  FOR  REQUIRED  MAINTENANCE  ACTIONS  BY  WUC 
C  NCSRA  —  CREW  SIZES  FOR  REQUIRED  MAINTENANCE  ACTIONS  BY  WUC 
COMMON/UCOM5/  RMTBF(74) 

C  RMTBF  —  MEAN  TIME  BETWEEN  FAILURES  BY  WUC 
COMMON/UCOM6/  RSTATS(4,20) 

C  RSTATS  --  STATISTICS  ON  CREW  USE  FOR  UNSCHEDULED  MAINTENANCE 
COMMON/UCOM7/  WAITING(82) 

C  WAITING  —  INDICATES  WHETHER  AN  AIRCRAFT  IS  WAITING  FOR  A  REPAIR  CREW 
COMMON/UCOM8/  RTRA(74,2) 

C  RTRA  --  REPAIR  TIME  FOR  A  REQUIRED  MAINTENANCE  ACTION  BY  WUC 
COMMON/UCOM9/  SORTAR(82,2) 

C  SORTAR  —  ARRAY  USED  TO  ORDER  AIRCRAFT  IN  EVENT  3 
COMMON/UCOMIO/  NCRWRQ( 74,13) ,NCRWSZ( 74 , 1 3 ) 

C  NCRWRQ  --  CREW  CODES  FOR  POSSIBLE  MAINTENANCE  ACTIONS  BY  WUC 
C  NCRWSZ  --  CREW  SIZES  FOR  POSSIBLE  MAINTENANCE  ACTIONS  BY  WUC 
COMMON/UCOMl 1 /  CRWPRB( 74 ,13), REPTIM( 74,13) 

C  CRWPRB  --  PROBABILITY  THAT  A  CREW  WILL  BE  USED  TO  REPAIR  A  POSSIBLE 
C  MAINTENANCE  ACTION 

C  REPTIM  “  REPAIR  TIME  FOR  A  POSSIBLE  MAINTENANCE  ACTION  BY  WUC 


oono  onorsnon 


C0MM0N/UC0M12/  IDTCC(32 , 7A ) , IDTCS(82 , 74 ) , IFAIL(82  ,74) 

C  IDTCC  “  CREW  CODES  OF  CREWS  DETERMINED  TO  REPAIR  CURRENT  FAILURES 
C  TO  AN  AIRCRAFT 

C  IDTCS  --  CREW  SIZES  OF  THE  CREWS  NEEDED  TO  REPAIR  CURRENT  FAILURES 
C  TO  AN  AIRCRAFT 

C  IFAIL  --  LIST  OF  CURRENT  FAILURES  TO  AN  AIRCRAFT 
COMMON/UCOM13/  DTRT(82,74) 

C  DTRT  —  REPAIR  TIMES  FOR  CURRENT  FAILURES  TO  AN  AIRCRAFT 
C  ARRAYS  ARE  DIMENSIONED  AS  FOLLOWS 
C  (74)  -  THERE  ARE  74  WUCS 

C  (74,2)  -  74  WUCS  BY  NUMBER  OF  REQUIRED  MAINTENANCE  ACTIONS  (THE  MODEL 

C  ONLY  USES  ONE  RMA  BUT  CAN  BE  EXPANDED  TO  MORE  THAN  ONE.) 

C  (74,13)  -  74  WUCS  BY  :IAXIMUM  OF  13  POSSIBLE  MAINTENANCE  ACTIONS 
C  (20,3)  -  THERE  ARE  20  TYPES  OF  SPECIALISTS  AND  3  SHIFTS 
C  (4,20)  -  4  TYPES  OF  STATISTICS  BY  20  SPECIALISTS 

C  STATISTICS 

1  -  NUMBER  OF  AIRCRAFT  CURRENTLY  WAITING  FOR  EACH  SPECIALIST 

2  -  TOTAL  NUMBER  OF  WAITS  FOR  EACH  SPECIALIST 

3  -  AVERAGE  AMOUNT  OF  TIME  AIRCRAFT  WAIT  FOR  EACH  SPECIALIST 

4  -  TOTAL  NUMBER  OF  USES  (IN  UNSCHEDULED  MAINTENANCE  ONLY)  FOR 
EACH  SPECIALIST 

(82,2)  -  THERE  ARE  82  AIRCRAFT  BY  2  ATTRIBUTES  USED  WHEN  SORTING 

(82,74)  -  82  AIRCRAFT  BY  74  POSSIBLE  FAILURES 
NCRDR  =  5 
NPRNT  =  6 
NTAPE  =  7 
NNSET  =  10000 
NWUC  =  74 
NCREWS  =  20 
CALL  SLAM 
STOP 
END 

ieieiiificific^'kitirk'kicieifiticicicificicif'tiieieieiciciiicicieifieic'kicicic'k'k'k'kicificikic'k’kic’k-k'kic'kic-kie-kic-k'kic'ic'k'k'k 

* 

SUBROUTINE  INTLC  * 

* 

C*********************************************************************** 

C  THIS  SUBROUTINE  READS  DATA  INTO  THE  ISHIFTS ,NRMA ,NPMA ,RMTBF ,RTRA, 

C  NCSRA,NCCRA,CRWPRB,REPTIM,NCRWSZ,  AND  NCRWRQ  ARRAYS.  IT  MODIFIES  THE 
C  DATA  FOR  ISHIFTS  AND  CRWPRB  TO  THE  FORM  NEEDED  IN  THE  MODEL.  IT  ALSO 
C  CREATES  THE  INITIAL  AMOUNT  OF  CREW  RESOURCES  IN  THE  SLAM  NETWORK. 
SUBROUTINE  INTLC 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET( 1 ) ) 

COMMON/SCOMl/  ATRIB( 100) ,DD( 100) ,DDL( 100) ,DTNOW, II ,MFA ,MSTOP ,NCLNR 
1 , NCRDR, NPRNT, NNRUN, NNSET, NTAPE, SS( 100 ),SSL( 100), TNEXT,TNOW,XX( 100) 


COMMON/UCOMl/  ITOTFL .NCREWS ,NWUC 
COMMON/UCOM2/  ISHFTS(20,3) 

COMMON/ UCOH3/  NFMAT(74),NPMA(74),NRMA(74) 

C0MM0N/UC0M4/  NCCRA(74,2) ,NCSRA{74,2) 

COMMON/UCOM5/  RMTBF(74) 

COMMON/UCOM8/  RTRA(74,2) 

COMMON/UCOMIO/  NCRWRQ( 74 ,13 ) ,NCRWSZ(74 , 13 ) 

COMMON/UCOMl 1/  CRWPRB( 74 , 13) ,REPTIM(74,13) 

DOUBLE  PRECISION  DCPACC 

OPEN(UNIT=NTAPE,FILE='MAINT. DAT', STATUS  =  'OLD') 
REWIND(UNIT=NTAPE) 

DO  100  I  =  1 ,NWUC 

READ(NTAPE,10)  RMTBF(l) ,NRMA( I) ,NPMA( l) 
READ(NTAPE,10) 

DO  150  J=1,NRMA(I) 

READ( NTAPE ,15)  RTRA( I , J ) ,NCSRA( I , J ) ,NCCRA ( I , J ) 
150  CONTINUE 

IF  (NRMA(I).GT.O)  READ(NTAPE,15) 

DCPACC  =0.0 

DO  200  J  =  1,NPMA(I) 

READ( NTAPE , 20 )  CRWPRB( I , J) ,REPTIM( I , J) , 

+  NCRWSZ(I,J) ,NCRWRQ(I,J) 

DCPACC  =  DCPACC  +  CRWPRB(I,J) 

CRWPRB(I,J)  =  DCPACC 
200  CONTINUE 

CRWPRB(I,NPMA(I))  =  1.0 

READ( NTAPE ,10 ,END=300) 

100  CONTINUE 

300  CLOSE(UNIT=NTAPE) 

OPEN(UNIT=NTAPE,FILE='CREW. DAT', STATUS  =  'OLD') 
REWIND(UNIT=NTAPE) 

DO  400  I  =  1,NCREWS 

READ( NTAPE, 30, END=401)  (ISHFTS(I,J) ,J=1 ,3) 

400  CONTINUE 

401  aOSE(UNIT=NTAPE) 

DO  500  I  =  1,NCREWS 

CALL  ALTER(I,ISHFTS(I,1)) 

500  CONTINUE 
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DO  600  I  =  1 ,NCREWS 

IMID  =  ISHFTS(I,1)  -  ISHFTS(I,3) 
IDAY  =  ISHFTS(I,2)  -  ISHFTS(I,1) 
ISWING  =  ISHFTS(I,3)  -  ISHFTS(I,2) 
ISHFTSd,!)  =  IMID 
ISHFTS(I,2)  =  IDAY 

ISHFTSd, 3)  =  ISWING 
600  CONTINUE 

10  F0RMAT(4X,F11.4,I4,I4) 

15  FORMATd7X,F10.6,6X,I2,I3) 

20  FORMATdlX,F6.3  ,F10.6,6X,I2,I3) 

30  FORMAT(5X,I2,I3,I3) 

RETURN 

END 


Cirk*****ic-k*******************'irk*1rk*-k1rk*********-)c**irk****irk*ie****-kie*ie**'k* 
SUBROUTINE  EVENT  (J) 

GO  TO  (1,2, 3, 4, 5, 6, 7, 8, 9),  J 

1  CALL  FAIL 
RETURN 

2  CALL  FREECREWS 

RETURN 

3  CALL  SCHEDULE 

RETURN 

4  CALL  DETCREWS 

RETURN 

5  CALL  TAKEOFF 

RETURN 

6  CALL  DAYSHIFT 

RETURN 

7  CALL  SWINGSHIFT 

RETURN 

8  CALL  MIDSHIFT 

RETURN 

9  CALL  TOPREFLIGHT 

RETURN 

END 
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c*********************************************************************** 

c  * 

SUBROUTINE  FAIL  * 

* 

(]*********★******■*****★********★*★************************************** 
C  THIS  SUBROUTINE  DETERMINES  WHETHER  A  FAILURE  HAS  OCCURRED  TO  ANY  OF 
C  THE  AIRCRAFT  SUBSYSTEMS  BASED  ON  A  PROBABILITY  OF  FAILURE  PER  FLIGHT. 

C  0.5  IS  ADDED  TO  THE  PROBABILITY  TO  USE  THE  MIDDLE  OF  THE  DISTRIBUTION 
C  GENERATING  RANDOM  NUMBERS.  THIS  IS  DONE  TO  INCREASE  THE  ACCURACY  OF 
C  THE  FAILURE  GENERATOR.  DOUBLE  PRECISION  IS  ALSO  USED  FOR  THIS  REASON. 
C  THIS  SUBROUTINE  ALSO  COUNTS  THE  NUMBER  OF  FAILURES  FOR  THIS  AIRCRAFT 
C  (NFAIL),  TOTAL  NUMBER  OF  FAILURES  (ITOTFL)  AND  BY  IRJC  (NFMAT). 
SUBROUTINE  FAIL 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,OSET( 1) ) 

COMMON/SCOMl/  ATRIB( 100 ) ,DD( 100) ,DDL( 100) ,DTNOW, II ,MFA,MSTOP,NaNR 
1 ,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOMl/  ITOTFL ,N CREWS ,NWUC 
COMMON/ U COM3 /  NFMAT ( 7 4 ) , NPMA ( 7 4 ) , NRMA (74) 

COMMON/UCOM5/  RMTBF(74) 

COMMON/UCOMl 2/  IDTCC(82 , 74) , IDTCS(82 , 74) , IFAIL( 82 , 74) 

DOUBLE  PRECISION  RNDNM.FLPCNT 

NPLANE  =  ATRIB(2) 

NFAIL  =  0 

DO  5  I  =  l.NWUC 

RNDNM  =  DRAND(9) 

FLPCNT=1.3Df00/RMTBF(I)+0.5D4-00 
IF  ((RNDNM. LT. 0.5). OR. (RNDNM. GT.FLPCNT))  GO  TO  5 
NFAIL  =  NFAIL  +  1 
NFMAT(I)  =  NFMAT(I)  +  1 
IFAIL(NPLANE, NFAIL)  =  I 
ITOTFL  =  ITOTFL  +1 
5  CONTINUE 

ATRIB(4)  *  NFAIL 


RETURN 


rjooo  ooooo 


THIS  SUBROUTINE  FREES  THE  RESOURCES  THAT  HAVE  BEEN  USED  FOR  A  REPAIR 


SUBROUTINE  FREECREWS 


COMMON/SCOMl/  ATRIB( 100 ) ,DD( 100) ,DDL( 100 ) ,DTNOW, II ,MFA ,MSTOP,NCLNR 
1 , N  CRDR , NPRNT ,NNRUN , NN  SET , NTAPE ,SS(100),SSL(100),TNEXT, TNOW , XX ( 1 00 ) 


J  =  ATRIB(6) 
K  =  ATRIB(7) 


CALL  FREE(J,K) 


ATRIB(5)=0.0 

ATRIB(6)=0.0 

ATRIB(7)=0.0 


RETURN 

END 


* 

SUBROUTINE  SCHEDULE  * 


C<r*****<f***'jWr  **********************  **********♦******♦♦*******<>**♦♦***♦** 

C  THIS  SUBROUTINE  SUBROUTINE  DOES  2  DIFFERENT  ACTIONS.  FIRST,  IT  SETS 
C  XX(13)  EQUAL  TO  THE  AVERAGE  VALUE  OF  FMC  AIRCRAFT  AS  COLLECTED  BY 
C  XX(1)  IN  THE  SLAM  NETWORK.  SECOND,  THIS  SUBROUTINE  REMOVES  ALL 
C  AIRCRAFT  IN  A  FILE  WAITING  FOR  TAKEOFF  AND  PLACES  THEM  INTO  THE  FILE 
C  FOR  SCHEDULING.  IT  THEN  STORES  IN  ARRAY  SORTAR  THE  ORDER  (BY  NUMBER 
C  OF  FLIGHT  HOURS)  OF  THE  AIRCRAFT  IN  THE  SCHEDULING  FILE.  FINALLY, 

C  THIS  SUBROUTINE  CALLS  THE  EVENT  TOPREFLIGHT. 

SUBROUTINE  SCHEDULE 
DIMENSION  NSET( 10000) 

COMMON  Q SET (10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET( 1 ) ) 

COMMON/SCOMl /  ATRI B( 100 ) , DD( 100 ) , DDL( 1 00 ) , DTNOW , I I , MFA , MSTOP , NCLNR 
1 ,NCRDR, NPRNT, NNRUN.NNSET, NTAPE, SS(100),SSL(100),TNEXT, TNOW, XX(IOO) 
COMMON /UCOM9/  SORTAR (8 2, 2) 


XX(13)  =  TTAVG(l) 


INDEX  »  0 

DO  10  I-  1,NNQ(7) 

CALL  RMOVEd  ,7,ATRIB) 
INDEX  -  INDEX  +  1 
SORTAR( INDEX , 1 ) -ATRIB( 2 ) 
SORTAR ( INDEX , 2 ) -ATR I B ( 3 ) 
CALL  FILEMdO.ATRIB) 
CONTINUE 

DO  20  1=  1,NNQ(3) 

CALL  RMOVEd, 3, ATRIB) 
INDEX  =  INDEX  +  1 
S0RTAR(INDEX,1)-ATRIB(2) 
SORTAR ( INDEX , 2 ) »ATR I B ( 3 ) 
CALL  FILEMdO.ATRIB) 
CONTINUE 

DO  30  1=  l.NNQdO) 

CALL  COPYd,  10, ATRIB) 
INDEX  =  INDEX  +  1 
SORTAR( INDEX , 1 )=ATRIB( 2) 
SORTAR( INDEX , 2 ) »ATRI B( 3 ) 
CONTINUE 


SWITCH  =0.0 
DO  50  1=2, INDEX 
IM1*I-1 

IF  (S0RTAR(IM1,2).GT.S0RTAR(I,2))  THEN 
T1=S0RTAR(IM1,1) 

T2=S0RTAR(IM1,2) 

SORTAR( IMl , 1 ) =SORTAR( 1,1) 

SORTAR ( IMl , 2 ) =SORTAR (1,2) 

SORTAR ( I, 1)=T1 
SORTAR(I,2)=T2 
SWITCH  =  1.0 
ENDIF 
CONTINUE 

IF  (SWITCH. GT. 0.0)  GO  TO  40 
CALL  SCHDL(9,0, ATRIB) 


RETURN 


o  o  o 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


********************************************** 


* 

* 

* 


SUBROUTINE  TO  PRE  FLIGHT 


* 

* 

* 


THIS  SUBROUTINE  MOVES  AIRCRAFT  FROM  THE  SCHEDULING  FILE  TO  THE 
QUEUE  FOR  PREFLIGHTS.  THE  AIRCRAFT  ARE  MOVED  IN  ORDER  OF  FLIGHT 
HOURS,  THE  HIGHEST  IS  MOVED  FIRST. 

CURRENTLY,  ALL  AIRCRAFT  ARE  ALLOWED  TO  FLY,  THUS  ALL  OF  THE  AIRCRAFT 
ARE  MOVED  AND  THE  CALL  OPEN(3)  STATEMENT  IS  USED  TO  OPEN  THE 
SCHEDULING  GATE.  IF  ALL  AIRCRAFT  ARE  NOT  ALLOWED  TO  FLY  THE  CODE  TO 
SELECT  THE  APPROPRIATE  AIRCRAFT  SHOULD  BE  INSERTED  AND  THE  CALL  OPEN 
STATEMENT  REMOVED. 

SUBROUTINE  TOPREFLIGHT 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET(1 ) ,QSET( 1)) 

COMMON/SCOMl/  ATRIB( 100) ,DD( 100) ,DDL( 100) ,DTNOW, II,MFA,MSTOP,NCLNR 
1 ,NCRDR ,NPRNT ,NNRUN ,NNSET ,NTAPE , SS( 100 ) , SSL( 100 ) , TNEXT , TNOW , XX( 100 ) 
COMMON/UCOM9/  SORTAR(82,2) 


DO  60  I  =  NNQ(10),1,-1 

NRANK  =  NFIND(l,10,2,0,SORTAR(I,l),.l) 
IF  (NRANK. GT.O)  THEN 

CALL  RMOVE ( NRANK, 10,ATRIB) 

CALL  FILEM(1,ATRIB) 

ENDIF 

60  CONTINUE 


CALL  OPEN(3) 


RETURN 

END 


C  THIS  SUBROUTINE  DETERMINES  WHICH  OF  SEVERAL  POSSIBLE  REPAIR  CREWS  WILL 


C  REPAIR  A  FAILURE. 

SUBROUTINE  DETCREWS 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET( 1 ) ) 

COMMON/SCOMl/  ATRIB( 100) ,DD( 100) ,DDL(100) ,DTNOW,II,MFA,MSTOP,NCLNR 
1, NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS( 100), SSL( 100), TNEXT, TNOW, XX( 100) 
COMMON/UCOMl/  ITOTFL ,NCREWS ,NWUC 
COMMON/UCOM3/  NFMAT(74) ,NPMA(74) ,NRMA(74) 

COMMON/UCOMIO/  NCRWRQ(74,13) ,NCRWSZ(74,13) 

COMMON/UCOMl 1 /  CRWPRB( 74,13) ,REPTIM( 74,13) 

COMMON/UCOMl 2/  IDTCC(82,74) ,IDTCS(82,74) ,IFAIL(82,74) 

COMMON/UCOMl 3/  DTRT(82,74) 
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NPLANE  =  ATRIB(2) 
NFAIL  =  ATRIB(4) 


DO  10  1  =  1,  NFAIL 
RNDNM  =  DRANDO) 

DO  20  J  =  1,NPMA(IFAIL(NPLANE,I)) 

JMl  =  J-1 

IF  ( ( RNDNM. LE . CRWPRB( IFAIL( NPLANE , I) , J ) ) . AND. 

+  (RNDNM.GT.CRWPRB(IFAIL(NPLANE,I),JM1)))  THEN 

IDTCC( NPLANE, I)  =  NCRWRQ( IFAIL(NPLANE,I) ,J) 
IDTCS(NPLANE,I)  =  NCRWSZ( IFAIL( NPLANE, I) ,J) 

DTRT( NPLANE, I)  =  REPTIM( IFAIL(NPLANE,I) ,J) 

ENDIF 

20  CONTINUE 

10  CONTINUE 

RETURN 

END 

C***»*************^Hr*<rtH^****************************ilr***<nllr********<riHr<>r** 

C  * 

C  SUBROUTINE  TAKEOFF  * 

C  * 

Cirk*1rk*'irkit*1e*1rk*it1t*1rk-)rkic1c*1t1c1rk1c1c**1titie1rk1rk1rk1rk1rk1tk1c****1rk*ic-k****1i1i*-k-k*-k'llric 

C  THIS  SUBROUTINE  MOVES  ONE  AIRCRAFT  TO  A  FILE  WHERE  THE  AIRCRAFT  CAN 
C  BEGIN  THE  FLIGHT  SEQUENCE. 

SUBROUTINE  TAKEOFF 

COMMON/SCOMl/  ATRIB(IOO)  ,DD(  100)  ,DDL(  100)  ,DTNOW, II,MFA,MSTOP,NaNR 
1 ,NCRDR ,NPRNT,NNRUN ,NNSET,NTAPE ,SS( 100) , SSL ( 100) ,TNEXT,TNOW,XX( 100) 

IF  (NNQ(7).GT.0)  THEN 

CALL  RMOVEd ,7,ATRIB) 

CALL  FILEM(8,ATRIB) 

ENDIF 


RETURN 

END 


C********************** AAA**  ft  k k-klcieifklrkirkirkitlrkiik  tthiflrirkle-klrklelclclelrle 


C  THIS  EVENT  ALTERS  THE  CREW  RESOURCES  FROM  THE  MIDNIGHT  SHIFT  TO  THE 
C  DAY  SHIFT  LEVELS. 

SUBROUTINE  DAYSHIFT 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET(1 )) 

COMMON/SCOMl/  ATRIB(  100)  ,DD(  100)  ,DDL(  100)  ,DTNOW,  II  ,MFA,MSTOP,NaNR 
1 ,NCRDR , NPRNT , NNRUN , NNSET , NTAPE , SS( 1 00 ) , SSL ( 1 00 ) , TNEXT , TNOW , XX( 100 ) 
COMMON/UCOMl/  ITOTFL ,NCREWS ,NWUC 
COMMON/UCOM2/  ISHFTS(20,3) 
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THIS  EVENT  ALTERS  THE  CREW  RESOURCES  FROM  THE  DAY  SHIFT  TO  THE  SWING 
SHIFT  LEVELS. 

SUBROUTINE  SWINGSHIFT 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

COMMON/SCOMl/  ATRIB(IOO) ,DD(100) ,DDL( 100) ,DTNOW, II,MFA,MSTOP,NCLNR 
1,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOMl/  ITOTFL,NCREWS,NWUC 
COMMON/UCOM2/  ISHFTS( 20 ,3 ) 

DO  100  I  =  1,NCREWS 

CALL  ALTER(I,ISHFTS(I,3)) 

)0  CONTINUE 

RETURN 

END 


SUBROUTINE  MIDSHIFT 


C  THIS  EVENT  ALTERS  THE  CREW  RESOURCES  FROM  THE  SWING  SHIFT  TO  THE 
C  MIDNIGHT  SHIFT  LEVELS. 

SUBROUTINE  MIDSHIFT 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

COMMON/SCOMl/  ATRIB(  100)  ,DD(  100)  ,DDL(  100)  ,DTNOW, II ,MFA,MSTOP,NaNR 
1,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOMl/  ITOTFL ,NCREWS ,NWUC 
COMMON/UCOM2/  1SHFTS(20,3) 

DO  100  I  =  1,NCREWS 

CALL  ALTER(I,ISHFTS(I,1)) 

100  CONTINUE 


RETURN 


111 
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* 

SUBROUTINE  ALLOC  * 

* 

C  THIS  SUBROUTINE  SEIZES  CREW  RESOURCES  (WHEN  AVAILABLE)  TO  FIX  AIRCRAFT 
C  IF  NO  RESOURCES  ARE  AVAILABLE  AND  THERE  ARE  MORE  AIRCRAFT  WAITING  FOR 
C  CREWS,  THE  AIRCRAFT  LAST  IN  THE  QUEUE  IS  REMOVED  AND  REPLACED  IN  THE 
C  QUEUE  TO  BE  CHECKED  FOR  CREW  AVAILABILITY.  THIS  OCCURS  UNTIL  ALL 
C  AIRCRAFT  HAVE  BEEN  CHECKED. 

SU  BROUTINE  ALLOC ( I , IFLAG ) 

DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

COMMON/SCOMl/  ATRIB(  100 )  ,DD(  100)  ,DDL(  100 )  ,DTNOW,  II  ,MFA  ,MSTOP  ,NaNR 
1,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 

IFLAG  -  0 

IF  (ATRIB(IO).GT.O.O)  THEN 

CALL  RMAWIP(NEED,ISIZE,RTIME) 

ELSE 

CALL  NOWIP(NEED,ISIZE,RTIME) 

ENDIF 

XX(ll)  =  XX(ll)  +  1.0 

IF  (ISIZE.EQ.O)  THEN 

IF  (XX(11).GE.NNQ(5))  THEN 
XX(ll)  =  0.0 
ELSE 

CALL  RMOVE(NNQ(5),5,ATRIB) 

CALL  FILEM(5,ATRIB) 

ENDIF 

ELSE 

CALL  SEIZE(NEED,ISIZE) 

ATRIB(5)=RTIME 
ATRIB(6)»NEED 
ATRIB(7)=ISIZE 
IFLAG  »  1 
ENDIF 

RETURN 

END 
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* 

* 

* 


SUBROUTINE  NOWIP 


* 

* 

* 


WHEN  CREWS  ARE  AVAILABLE  THIS  SUBROUTINE  PROVIDES  ALLOC  WITH  THE  CREW, 
ITS  SIZE,  AND  THE  REPAIR  TIME  NEEDED  TO  FIX  A  FAILURE. 

WHEN  CREWS  ARE  NOT  AVAILABLE  THE  SIZE  VARIABLE  IS  A  ZERO. 

THIS  SUBROUTINE  GETS  THIS  INFORMATION  FROM  AVAIL. 

IF  THIS  REPAIR  ACTION  WILL  FINISH  FIXING  A  FAILURE  THIS  SUBROUTINE 
CAUSES  THE  ARRAYS  THAT  STORE  FAILURE  DATA  TO  BE  REDUCED. 

SUBROUTINE  NOWIP(NEED, ISIZE,RTIME) 

DIMENSION  NSET(IOOOO) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET(1 ) ) 

COMMON/SCOMl/  ATRIB(IOO) ,DD( 100) ,DDL( 100) ,DTNOW, II,MFA,MSTOP,NaNR 
1 ,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOM3/  NFMAT(74) ,NPMA(74) ,NRMA(74) 

COMMON/UCOM12/  IDTCC(82 ,74) ,IDTCS(82 ,74) ,IFAIL(82 ,74) 


NEED  =  0 
ISIZE  =  0 
RTIME  =0.0 


NPLANE  =  ATRIB(2) 

NFAIL  =  ATRIB(4) 

DO  10  1  =  1,  NFAIL 

MA  =  IFAIL( NPLANE, I) 

IF  (NRMA(MA).GT.O)  THEN 
NACT  =  0 
NEXTRA  =  1 
ELSE 

NACT  =  1 
NEXTRA  »  0 
ENDIF 


CALL  AVAIL( NPLANE , NACT, I ,MA , NEXTRA , ICC , ICSZND ,REPLEN) 
IF  (ICSZND. GT.O)  GOTO  30 

10  CONTINUE 

CALL  CREWSTATS ( NPLANE , NFAIL , ICC ,0 ) 

RETURN 


30  CALL  CREWSTATS (NPLANE, NFAIL, ICC, ICSZND) 


IF  (NACT.EQ.l)  CALL  REDUCE(l) 

NEED  =  ICC 
ISIZE  *  ICSZND 
RTIME  =•  REPLEN 


RETURN 

END 


C  *  * 

C  *  SUBROUTINE  RMAWIP  * 

C  *  * 

C  WHEN  CREWS  ARE  AVAILABLE  THIS  SUBROUTINE  PROVIDES  ALLOC  WITH  THE  CREW, 
C  ITS  SIZE,  AND  THE  REPAIR  TIME  NEEDED  TO  FIX  A  FAILURE. 

C  WHEN  CREWS  ARE  NOT  AVAILABLE  THE  SIZE  VARIABLE  IS  A  ZERO. 

C  THIS  SUBROUTINE  GETS  THIS  INFORMATION  FROM  AVAIL. 

C  IF  THIS  REPAIR  ACTION  WILL  FINISH  FIXING  A  FAILURE  THIS  SUBROUTINE 
C  CAUSES  THE  ARRAYS  THAT  STORE  FAILURE  DATA  TO  BE  REDUCED. 

SUBROUTINE  RMAWIP(NEED , ISIZE , RTIME) 

DIMENSION  NSET(IOOOO) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET(l) ,QSET(1)) 

COMMON/SCOMl/  ATRIB(  100)  ,DD(100)  ,DDL(  100)  ,DTNOW,II,MFA,MSTOP,NaNR 
1 ,NCRDR,NPRNT,MNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOM3/  NFMAT(7A) ,NPMA(7A) ,NRMA(74) 

NEED  =  0 
ISIZE  =  0 
RTIME  =  0 

NPLANE  =  ATRIB(2) 

MA  =  ATRIB(9) 

IFLNUM  =  ATRIB(ll) 

NEXTRA  =  ATRIB(IO)  +  1 

IF  ( NEXTRA. GT.NRMA(MA))  THEN 
NACT  =  1 
ELSE 

NACT  =  0 
ENDIF 

CALL  AVAIL ( NPLANE , NACT , IFLNUM, MA , NEXTRA , ICC , ICSZND , REPLEN) 

CALL  CREWSTATS(NPLANE,1 , ICC, ICSZND) 


I 
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IF  ((NACT.EQ.l).AND.(ICSZND.GT.O))  CALL  REDUCE( IFLNUM) 


NEED  =  ICC 
ISIZE  =  ICSZND 
RTIME  =  REPLEN 

RETURN 

END 


**************************************** 


SUBROUTINE  AVAIL 


C  THIS  SUBROUTINE  PROVIDES  NOWIP  AND  RMAWIP  WITH  THE  CREW  CODE,  ITS  SIZE 
C  AND  REPAIR  TIME  NEEDED  TO  FIX  A  FAILURE.  IF  NOT  ENOUGH  CREWS  ARE 
C  AVAILABLE  THE  CREW  SIZE  IS  SET  TO  ZERO. 

SUBROUTINE  AVAIL ( NPLANE ,NACT , IFLNUM, MA ,NEXTRA ,NEED , ISIZE , RTIME) 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET( 1) ) 

COMMON/SCOMl/  ATRIB( 100) ,DD( 100) ,DDL( 100 ) ,DTNOW, II ,MFA ,MSTOP ,NCLNR 
1 ,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COHMON/UCOM4/  NCCRA(74 ,2) ,NCSRA(74 ,2) 

COMMON/UCOM8/  RTRA(74,2) 

COMMON/UCOM12/  IDTCC(82,74) ,IDTCS(82,74) ,IFAIL(82,74) 
COMMON/UCOM13/  DTRT(32,74) 

IF  (NACT.GT.O)  THEN 

ICC  =  IDTCC(NPLANE, IFLNUM) 
lAVAIL  =  NNRSC(ICC) 

ICSZND  =  IDTCS( NPLANE, IFLNUM) 

REPLEN  =  DTRT( NPLANE, IFLNUM) 

ELSE 

ICC  =  NCCRA(MA,NEXTRA) 
lAVAIL  =  NNRSCUCC) 

ICSZND  =  NCSRA(MA,NEXTRA) 

REPLEN  =  RTRA(MA,NEXTRA) 

ENDIF 


IF  ( (lAVAIL. LT. ICSZND). OR. (lAVAIL. EQ.O))  THEN 
ICSZND  =  0 
ELSE 

IF  (NACT.EQ.l)  THEN 
ATRIB(9)  =0.0 
ATRIB(IO)  =  0.0 
ATRIE(ll)  =  0.0 
ELSE 

ATRIB(9)  =  MA 
ATRIB(IO)  =  NEXTRA 
ATRIB(ll)  =  IFLNUM 
ENDIF 
ENDIF 
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NEED  =  ICC 
ISIZE  =  ICSZND 
RTIME  =  REPLEN 

RETURN 

END 

C  **************************************** 

c  *  * 

c  *  SUBROUTINE  CREWSTATS  * 

C  *  * 

C  ****************************************** 

C  THIS  SUBROUTINE  COLLECTS  STATISTICS  ON  HOW  MANY  TIMES  CREWS  ARE  USED 
C  FOR  UNSCHEDULED  REPAIRS.  IT  ALSO  COLLECTS  HOW  MANY  TIMES  CREWS  WERE 
C  NEEDED  BUT  NOT  AVAILABLE  AND  THE  AVERAGE  WAITING  TIME  FOR  THESE  CREWS 
SUBROUTINE  CREWSTATS(NPLANE,NCHECK, NEED, ISIZE) 

DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET( 1 ) ) 

COMMON/SCOMl/  ATRIB(IOO) ,DD(100) ,DDL(100) ,DTNOW,II,MFA,MSTOP,NCLNR 
1 ,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100),SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOM3/  NFMAT(74) ,NPMA(74) ,NRMA(74) 

C0MM0N/UC0M4/  NCCRA( 74 , 2 ) ,NCSRA( 74 , 2 ) 

COMMON/UCOM6/  RSTATS(4,20) 

COMMON/UCOM7/  WAITING(82) 

COMMON/UCOMl 2/  IDTCC(82,74), IDTCS (82, 74), I FAIL (82, 74) 

DO  100  I  «  1,NCHECK 

IF  (NCHEClC.EQ.l)  THEN 
ICC  =  NEED 
ELSE 

MA  =  IFAIL(NPLANE,I) 

IF  (NRNA(MA).GT.O)  THEN 
ICC  =  NCCRA(MA,1) 

ELSE 

ICC  =  IDTCC(NPLANE,I) 

ENDIF 

ENDIF 

IF  (ISIZE. GT.O)  THEN 

IF  (WAITING(NPLANE).GT.O.O)  THEN 

RSTATSd ,ICC)  =  RSTATS(1,ICC)  -  1 
RSTATS(3,ICC)  =  RSTATS(3 ,ICC)  -  TNOW 
ENDIF 
ELSE 

IF  (WAITING(NPLANE).EQ.O.O)  THEN 

RSTATSd, ICC)  “  RSTATSd, ICC)  +  1.0 
RSTATS(2,ICC)  =  RSTATS(2,ICC)  +  1.0 
RSTATS(3,ICC)  =  RSTATS(3,ICC)  +  TNOW 
ENDIF 
ENDIF 

100  CONTINUE 


'.V 
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IF  (ISIZE.GT.O)  THEN 

WAITING (NPLANE)  =0.0 
RSTATS(4,NEED)  =  RSTATS(4 .NEED)  +  1.0 
ELSE 

WAITING(NPLANE)  =1.0 
ENDIF 

RETURN 

END 


C  ************************************ 

c  *  * 

C  *  SUBROUTINE  REDUCE  * 

C  *  * 

C  ************************************ 

C  THIS  SUBROUTINE  REDUCES  THE  ARRAYS  THAT  STORE  WHICH  AIRCRAFT  HAVE 
C  WHAT  FAILURES,  WHO,  HOW  MANY,  AND  HOW  LONG  IT  TAKES  TO  FIX  THEM. 

C  IT  ALSO  DECREMENTS  THE  ATTRIBUTE  WHICH  STORES  HOW  MANY  FAILURES 
C  THE  AIRCRAFT  HAS. 

SUBROUTINE  REDUCE( INDEX) 

DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  ( NSET( 1 ) ,QSET(1 ) ) 

COMMON/SCOMl/  ATRIB( 100) ,DD( 100) ,DDL( 100 ) ,DTNOW, II ,MFA ,MST0P ,NCLNR 
1 ,NCRDR ,NPRNT ,NNRUN ,NNSET ,NTAPE , SS(100 ) , SSL(100 ) ,TNEXT , TNOW , XX( 100 ) 
COMMON/UCOM12/  IDTCC(82,74) ,IDTCS(82,74) ,IFAIL(82,74) 
COMMON/UCOM13/  DTRT(82,74) 


100 


NPLANE  =  ATRIB(2) 

NFAIL  =  ATRIE(4) 

IF  (INDEX.lt. NFAIL)  THEN 
INDXPl  =  INDEX  +  1 
DO  100  I  =  INDXPl  .NFAIL 
IMl  =  I  -  1 
IFAIL(NPLANE,IM1)  = 
IDTCC(NPLANE,IM1)  = 
IDTCS( NPLANE, IMl)  = 
DTRT( NPLANE, IMl)  = 
CONTINUE 
ENDIF 

IFAIL(NPLANE, NFAIL)  =  0 
IDTCC( NPLANE, NFAIL)  =  0 
IDTCS( NPLANE, NFAIL)  =  0 
DTRT(NPLANE, NFAIL)  =0.0 


IFAIL< NPLANE, I) 
IDTCC( NPLANE, I) 
IDTCSOlPLANE.l) 
DTRT( NPLANE, I) 


ATRIB(4)  =  NFAIL  -  1 


RETURN 

END 
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no-MS?  121 

UNCLASSIFIED 


A  SIHULATION  NODEL  OF  THE  T-4SA  AIRCRAFT  FOR 
AVAILASILITV  AND  SORTIE  PROJECTIONS! U>  AIR  FORCE  INST 
OF  TECH  URilHT-PATTERSON  AF8  OH  R  A  FOLEV  ET  AL. 

DEC  09  AFIT/MR/ENS/8SD-E  F/G  3/9 


NL 


non 


- •  V  ^  I.' 


Qi:icieicicicicisifieifi(icicifie'kiti:iei:icicicici:ic'icicic'ic'k’k-kis'kic'ki:is-k’kicig'k'kic'k’ic'k'k'kic’k'kie'k'!c'k'kic'k'k'k’kicisificicie 

* 

SUBROUTINE  OTPUT  * 

* 


C  THIS  SUBROUTINE  PRINTS  OUT  USER  DEFINED  OUTPUT 
SUBROUTINE  OTPUT 
DIMENSION  NSET( 10000) 

COMMON  QSET( 10000) 

EQUIVALENCE  (NSET( 1 ) ,QSET( 1 ) ) 

COMMON/SCOMl/  ATRIB(IOO) ,DD( 100) ,DDL( 100) ,DTNOW, II ,MFA,MSTOP,NCLNR 
1 ,NCRDR,NPRNT,NNRUN,NNSET,NTAPE,SS(100) ,SSL(100),TNEXT,TNOW,XX(100) 
COMMON/UCOMl/  ITOTFL , NCREWS ,NWUC 
COMMON/UCOM3/  NFMAT( 74 ) ,NPMA( 74) ,NRMA( 74 ) 

COMMON/UCOM4/  NCCRA( 74 ,2 ) ,NCSRA( 74 ,2) 

COMMON/UCOM5/  RMTBF(74) 

COMMON/UCOM6/  RSTATS(4,20) 

COMMON/UCOM8/  RTRA(74,2) 

COMMON/UCOM9/  SORTAR(82,2) 

COMMON/UCOMIO/  NCRWRQ( 74 ,13) ,NCRWSZ( 74 ,13) 

COMMON/UCOMl 1 /  CRWPRB< 74 , 1 3 ) ,REPTIM( 74 , 1 3 ) 

COMMON/UCOMl 2/  IDTCC(82 ,74) , IDTCS(82 ,74) , IFAIL(82 ,74) 
COMMON/UCOM13/  DTRT(82,74) 


CHARACTER* 5  CREW(20) 

DATA  CREW/'R431  ','F431  ','P431 
+'S426  ','A426  ','1426  ','F426  ' 
+'S4231','S4232','S3280','S3281' 


','S4270','S4274','S4275','S4271', 
,'W431  ','S4233','S4230','S4234', 
,'S325  '/ 


DIMENSION  EXPT(74) 

DOUBLE  PRECISION  FLPCNT 

WRITE(NPRNT,100) 

100  FORMATdX, 'OBSERVED  FAILURES  BY  WUC') 
yRITE(NPRNT,101)  (NFMAT( I) , 1=1 ,NWUC) 

101  FORMAT(10(3X,I5)) 

CHISQ  =0.0 
SUMFLS  =0.0 
DO  10  I  =  1 ,NWUC 

FLPCNT=1.3D+00/RMTBF(I) 

EXPT(I)  =  XX(3)  *  FLPCNT 
SUMFLS  =  SUMFLS  +  EXPT(I) 

CHISQ  =  CHISQ  +  ((NFMAT(I)-EXPT(I))**2)/EXPT(I) 
10  CONTINUE 

WRITE(NPRNT,101) 

WRITE(NPRNT,120) 

120  FORMATdX, 'EXPECTED  NUMBER  OF  FAILURES  BY  WUC') 
WRITE( NPRNT ,121)  ( EXPT( I ) , 1=1 ,NWUC) 

121  FORMATdO(F7.3 ,1X)) 

WRITE( NPRNT, 101) 
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WRITE(NPRNT,130)  CHISO 

FORMATdX, 'CHI-SQUARE  SFATISriC  FOR  FAILURE  CEUER  A  FOR  • 

WRITE(NPRNT,101) 

l/RITE(NPRNT,140)  XX(3) 

FOR.XAK  IX, 'TOTAL  NUMBER  OF  FLIGHTS  IS', FIG. 1) 
URITE(NPRNT,150)  ITOTFL 

FORMAT( IX, 'OBSERVED  NUMBER  OF  FAILURES' , 16 ) 

V.’RITE(HPRNT,160)  SUNFLS 

FORMATdX, 'EXPECTED  NUMBER  OF  FAILURES' ,  F8 . 2  ) 

DO  20  I=1,NWUC 
EXPT(I)  =0.0 
NFMATd)  =  0 
CONTINUE 

DO  40  I  =  1  ,NCREWS 

IF  (RSTATSd  ,I).EQ.0.0)  GO  TO  40 
NCW  =  RSTATS( 1,1) 

DO  30  J=1 ,NCW 

RSTATS(3,I)  =  RSTAFS(3,I)  -  TNOW 
CONTINUE 
CONTINUE 

DO  50  I  =  1 ,NCREWS 

IF  (RSTATS(2,I).EQ.0.0)  GO  TO  50 
RSFATS(3,I)  =  -RSTATS(3,I)/RSTATS(2,I) 

CONTINUE 

WRITE(NPRNT,170) 

FORMATdX, 'CREWS') 

WRITE(NPRNT,130)  (CREWd),I=l  ,NCREWS) 

FORMATdO(3X,A5)) 

WRITE(NPRNT,190) 

FORMAFdX, 'RESOURCE  STATS:  NUMBER  GF  FII^'ES  USED  TO  RFPAIR' ) 
v;RITF,(NPRNT,200)  (  R3TAFS(4  ,  l) .  1  =  1  ,NCREWS) 

FORMATdO(F8.1)) 

WRirr(N?RNT,310) 

FORMAFdX/RESOURCE  STA  FS :  NUMBER  CURRENFLi’  WAITING') 
WRITE(NPRNT,200)  (RSTATSI 1 , I) , 1=1 .NCREWS) 

v;rite(nprnt,210) 

FORMATdX, 'RESOURCE  STAFS:  FOIAL  NUMBER  OF  WAITS') 
'..■RITE(NPRNT,200)  (RSFAFSTU  ,  I)  ,  1  =  1  ,NCRr.WS) 

WRrFE(NPRNT,220) 

FORMATdX, 'RESOl'RCE  STAFS:  AVERACIE  WAITING  IIHE') 
WRITE(NPRNT,230)  ( RSF A FS( 3 ,  I ) , 1  =  1 .SCREWS) 

FORMATdO(FS,  1 ) ) 


Table  B.l 


Work  Unit  Code  Data 


wuc 

Nomenclature 

MTBF 

MTTR 

hours 

hours 

1  lA 

fuselage,  forward  section 

468 .1669 

3.587277 

1  1  B 

fuselage,  center  section 

1168.2249 

0.591818 

lie 

wing  assembly 

1642.0352 

1.11821 

1  ID 

empennage 

810 .3720 

2.391503 

HE 

engine  nacelle 

734.2144 

1.788432 

IIF 

fuselage,  aft  section 

4290.0000 

0.933081 

12A 

cockpit 

722.2833 

3.913892 

12B 

canopy 

623.2855 

4.386283 

1  2C 

ejection  seat  system 

418.9456 

2.656802 

13A 

main  landing  gear 

113.0546 

2.219961 

13B 

nose  landing  gear 

234.0806 

1.614848 

13C 

brake  system 

912.4177 

2.401157 

1  3D 

landing  gear  control  system 

2164.4912 

2 .18305 

13E 

auxiliary  landing  gear  extension 

2164.5000 

1.528385 

1  3F 

nose  wheel  steering 

1237.0099 

2.90717 

13G 

landing  gear  warning 

2134.4690 

1.54069 

14A 

pilot  controls 

2840.9060 

2 .298255 

14B 

roll  control 

853.9690 

3.93062 

14C 

pitch  control 

977.0434 

4.374974 

1  4D 

yaw  control 

806.2630 

1.475129 

14E 

trailing  edge  system 

793 .1480 

2.655312 

14F 

speed  brake  system 

1117.3239 

2.416026 

23A 

engine  ( core ) 

359.2000 

2.712 

2  3N' 

1  gn  i  1 1  on/ e lec t r 1 ca 1  system 

1628.6646 

1 .332284 

23P 

engine  lubrication  system 

3246.7517 

1.681402 

2  30 

main  fuel  system 

1329.8004 

4.136834 

2  3R 

engine  instrumentation  system 

346 .3588 

3 . 534767 

23  S 

starting  system 

8271 .2920 

2.059589 

23T 

engine  control  system 

305.9016 

4 . 594358 

23  V 

built-up  engine 

2238.314/ 

2. 138305 

4  1  A 

cockpit  air  temp  control  system 

847.4592 

1 .884295 

4  1  B 

air  conditioning 

1519.7465 

3.095565 

•4  1  C 

pressurization 

1736.1019 

3 .778195 

4  1  D 

bleed  air  system 

671.2315 

2.992214 

4  1  E 

anti-ice  system 

1126.0149 

0.842065 

4  1  F 

windshield  de-ice  system 

2053.6931 

1.207701 

4  1  r; 

d  e  f  0  g 

4032 . 2629 

1 .65441 

4  1  >i 

ram  air 

1  5290 . 5068 

1 .84326 

■4  1  J 

avionic  equipment  cooling 

3581.6770 

0 .976202 

42A  primary  DC  power  system 

42B  AC  power  system 

42C  DC  emergency  system 

42D  external  power  system 

42E  AC/DC  distribution  system 

44A  exterior  lighting  system 

44B  interior  lighting  system 

44C  caution  advisory  system 

45A  hydraulic  power  generation  system 

45C  hydraulic  indicator  system 

46a  fuel  storage  Installation 

46B  fuel  vent  installation 

46C  fuel  quantity  indicating  system 

46D  fuel  feed  system 

46E  ground  refuel  system 

46F  fuel  precheck  and  management  system 

47A  LOX  supply  system 

49A  fire  detection  system 

51A  flight  instruments 

51B  navigation  instruments 

51C  HARS  ,AN/ASN-129 

51D  pitot-static  system 

51E  cockpit  pressure 

52A  stability  augmentation  system 

55A  air  data  record  system 

62A  VHF/AM  communications  system 

63A  UHF  communications  system 

64A  intercommunications  system 

65A  transponder  set,  AN/ APX- 100( V) 

71A  VOR/ILS/MB  system  ARN-127 

71B  TACAN  system  AN/ARN-118 

91A  pilots  emergency  equipment 

91B  crash  position  indicator  system 

97A  canopy  removal  system 

97B  ejection  seat  removal  system 

(compiled 


711.1420 

1.707457 

2233.6328 

1.352932 

750.0913 

0 .730963 

13422.8184 

1.469011 

4905.9395 

7.621387 

384.4512 

0.788592 

779.8432 

0.629052 

2500.0046 

0.787688 

589.4842 

2.311582 

3601.2700 

1.178073 

30303.0020 

12.662677 

24390.2441 

5.105008 

3579.1226 

4.498434 

2077.2627 

3.469134 

12048 .1924 

4 .152074 

6655.1367 

1.867408 

422.6771 

2.029429 

1973.1656 

2.68424 

109.9207 

2.141921 

230.9148 

2.285361 

590.8571 

5 .197813 

1436.3468 

7.002 

4672.8999 

0.735 

1218.9034 

2.215331 

710.0895 

0.691788 

286.4782 

2.130949 

276.0772 

1.971631 

681.1048 

2.153563 

590.5653 

1  .513261 

491.9169 

2.436006 

528.3670 

2.623565 

2418.4241 

1.640232 

18222.2124 

1  .  74 

16666.6621 

0  .605119 

4739.3354 

3.411 

from  8 : A- ] 

[  to  A-38) 

Table  B.2  is  the  file  called  "crew.dat”  used  in  the 
T-46A  model.  The  specialty  code  corresponds  to  the  code  used 
to  label  the  resources  in  the  SLAM  code.  The  SLAM  resource 
code  is  the  ATC  LCOM  crew  code  with  the  alpha  character 
placed  first  rather  than  fourth.  The  fifth  character  is 
retained  if  needed  to  uniquely  identify  the  specialist.  For 
example,  the  LCOM  crew  code  427SS  becomes  the  SLAM  resource 
code  S4275.  The  shift  sizes  are  in  the  order  of  midnight 
shift,  day  shift,  and  swing  shift. 


Table  B.2 


Work  Center  Data 


Specialty 

Shift 

size 

Work  centers 

Code 

M 

D 

S 

R431 

0 

4 

2 

T-37 

repair  and  reclamation 

F431 

3 

38 

18 

T-37 

f  Im 

P431 

5 

5 

0 

T-37 

inspection 

S4270 

0 

2 

2 

machine 

S4274 

2 

2 

2 

metals 

processing 

S4275 

2 

6 

6 

structural  repair 

S4271 

1 

1 

1 

corrosion  control 

S426 

3 

3 

3 

T-37 

jet  engine  shop 

A426 

3 

2 

3 

T-37/38 

accessory  repair 

T426 

3 

3 

0 

T-37/38 

test  cell 

F426 

3 

9 

9 

T-37 

flight  line  support  unit 

W431 

2 

4 

3 

T-37/38 

wheel  &  tire 

S4233 

4 

4 

2 

T-37/38 

fuel  systems 

S4230 

4 

6 

6 

T-37/38 

electrical  systems 

S4234 

6 

6 

6 

T-37/38 

pneudraulics  systems 

S4231 

4 

4 

4 

T-37/38 

environmental  systems 

S4232 

2 

2 

2 

T-37/38 

egress  systems 

S3280 

2 

2 

2 

T-37/38 

radio  and  radar  repair 

S328  1 

2 

4 

4 

T-37/38 

radio  and  radar  repair 

S325 

6 

10 

10 

T-37/38  auto  flight  control 

(compiled  from  3:3-1  to  3-27) 

Table  B.3  contains  a  list  of  the  scheduled  maintenance 


as  proposed  by  Fairchild  Republic.  It  also  contains  the 
suggested  frequency  and  an  estimate  of  the  maintenance 
manhours  Involved. 


Table  B.3 

Scheduled  Maintenance  Actions 

Task  Nomenclature  MTBSM^  MMH^ 


special  Inspection 

rsvr,  master  cyl  brk 

25 

0.266667 

spectro  analysis 

basic  engine  (F109) 

25 

0.375002 

ADR  data  extraction 

recorder,  airborne  data  39 

0.083333 

clean/vacuum  Interior 

cockpit 

60 

0.266667 

special  capacitance 

battery  assembly 

60 

0 . 533333 

ADR  data  extraction 

recorder,  airborne  data  66. 

7  0.083333 

special  Inspection 

regulator  oil  demand 

120 

2.112 

replace  oil  filter 

basic  engine  (F109) 

150 

0.58497 

washing  A/C 

al rf  rame 

180 

12.8 

lube  due  to  washing 

a  1 r f  rame 

180 

0.533333 

ground  handling  A/C 

200  special  Inspection 

200 

0.96 

phase 

300 

21.163869 

lubrication 

300 

4.349602 

Insp/repack  kit 

survival  kit  assy 

300 

4.544 

corrosion  prevention 

airframe 

360 

3.2 

phase 

600 

2.235466 

Insp/repack  chute 

parachute  system  L/R 

720 

8.0 

replace  brushes 

star  ter/ generator 

1000 

10.8224 

engine  HSI 

basic  engine  (F109) 

1200 

10.0 

engine  bulld-up 

basic  engine  (F109) 

1200 

4.8 

nondestructive  Insp 

alrf  rame 

1500 

48.0 

lubrication 

starter/generator 

2000 

10.8224 

1  Mean  time  between  scheduled  maintenance  In  flight  hours. 

2  Maintenance  manhours. 

(compiled  from  8:4-10  to  4-15) 


Table  B.A  Is  the  file  "maint.dat”  used  In  the  T-46A 
model.  There  are  two  types  of  records  in  this  file.  Type 
one  records  contain  a  WUC,  the  MTBF  for  that  WUC  and  then  the 
number  of  required  and  possible  maintenance  actions.  Type 
two  records  contain  a  WUC,  a  node  label,  a  probability  of 
that  maintenance  action  being  taken  (only  for  possible 
maintenance  actions),  the  MTTR,  specialty  code  required,  crew 
size  required  and  a  numeric  code  for  the  specialty  code. 
Figure  B.l  contains  examples  of  the  record  types  from  WUC 


Type  1  record 

WUC  MTBF  NRMA 

IID  810.372  1 

NPMA 

5 

Type  2  record 

WUC  Node 

Probabll 1 ty 

MTTR 

Specialty 

Crew  Size 

Crew 

Label 

of  Selection 

Code 

Required 

Code 

IID  llAMl 

0.196  2 

.391503 

427S0 

1 

4 

Figure  B.l.  Example  of  Record  Types 


The  numeric  code  is  dependent  on  the  order  in  which  all 
resources  are  declared  in  the  SLAM  code.  The  model  has  been 
set  up  with  the  crew  resources  listed  first.  Since  there  are 
20  specialty  codes  the  crew  codes  range  from  1  to  20.  If  the 


order  in  which  resources  are  declared  in  the  SLAM  code 


changes,  the  numeric  code  must  be  changed.  Changing  the 
order  also  Impacts  the  shift  change  events  In  the  FORTRAN 
coding  where  the  loop  counter  Is  used  as  the  SLAM  resource 
code  (l.e.  the  crew  code).  The  node  label  Is  the  node  label 
from  the  appropriate  LCOM  model.  Nodes  with  a  W  as  the  first 
letter  of  the  label  are  from  the  ASD  network.  The  specialty 
codes  In  this  table  correspond  to  the  specialty  codes  that 
were  taken  from  the  ATC  LCOM  Final  Report. 

The  records  are  arranged  In  sets  In  the  file.  Each  set 
of  records  represents  a  repair  node  In  an  LCOM  network.  The 
first  record  of  each  set  Is  a  type  1  record.  Adding  the 
number  of  required  and  possible  maintenance  actions  Indicates 
how  many  type  2  records  follow.  For  example,  examine  the 
node  associated  with  repair  of  WUC  IID.  The  first  of  the  set 
(underlined  on  page  96)  shows  there  Is  one  required 
maintenance  action  and  five  possible  maintenance  actions. 
Therefore,  there  are  six  type  2  records  In  the  set  for  node 
IID.  The  first  of  these  records  Indicates  that  three  431F7 
repair  personnel  are  needed  for  0.8  hours.  Next  are  the 
records  for  the  five  possible  actions  that  will  finish  the 
repair  of  WUC  IID.  Each  of  these  has  a  probability  of 
occurring  but  only  one  will  be  done.  Estimates  of  how  long 
different  crews  would  take  to  repair  the  failure  are  not 
available.  Thus,  the  Fairchild  Republic  estimate  of  the  MTTR 
Is  used  as  the  repair  time  for  each  crew.  Figure  B.2  shows  a 
graphic  representation  of  node  IID. 
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Figure  B.2.  Repair  Node  IID 


Table  3.4 

Unscheduled  Maintenance  Network 


1  lA 

468  . 

1669 

0  6 

llA 

1 1  5M1 

0 .148 

3.587277 

427S0 

1 

4 

llA 

1 15M2 

0.011 

3.537277 

427S1 

1 

7 

llA 

115M3 

0.538 

3.587277 

427S5 

1 

6 

llA 

1 15M4 

0.130 

3.587277 

431F7 

1 

2 

llA 

115R1 

0.111 

3.587277 

431F7 

1 

2 

llA 

1  15R2 

0.062 

3 . 587277 

431R7 

1 

1 

IIB 

1168  . 

2249 

0  6 

IIB 

1 15M1 

0.148 

0.591818 

427S0 

1 

4 

IIB 

1 15M2 

0.011 

0.591818 

427S1 

1 

7 

IIB 

1 15M3 

0 .538 

0.591318 

427S5 

1 

6 

IIB 

115M4 

0.130 

0.591818 

431F7 

1 

2 

115R1 

0.111 

0 . 591818 

431F7 

1 

2 

115R2 

0.062 

0.591813 

431R7 

1 

1 

1642  . 

0352 

0  6 

1 13R1 

0.018 

1.118210 

423S3 

2 

13 

1 13M1 

0.020 

1 .118210 

427S1 

1 

7 

113M2 

0.628 

1.118210 

427S5 

1 

6 

113M3 

0 .106 

1 .118210 

431F7 

1 

2 

113R2 

0.106 

1.118210 

431F7 

1 

2 

113M4 

0.122 

1.118210 

427S0 

1 

4 

810  . 

3720 

1  5 

114D 

0.8 

431F7 

3 

2 

114M1 

0 .196 

2.391503 

427S0 

1 

4 

1 14M2 

0.067 

2.391503 

427S1 

1 

7 

114M3 

0.512 

2.391503 

427S5 

1 

6 

1 14M4 

0 .135 

2.391503 

431F7 

1 

2 

114R1 

0.090 

2.391503 

431F7 

1 

2 

734  . 

2144 

0  6 

113R1 

0.018 

1  .  738432 

423S3 

2 

13 

1 13M1 

0.020 

1.788432 

427S1 

1 

7 

113M2 

0.628 

1.733432 

427S5 

1 

6 

1  13M3 

0 .106 

1  .  738432 

431F7 

1 

2 

1 1  3  R  2 

0 . 106 

1  .  788432 

431F7 

1 

2 

1  13M4 

0.122 

1.788432 

427S0 

1 

4 

4290  . 

0000 

0  6 

llj’Al 

0.148 

0.933081 

427S0 

1 

4 

1 15M2 

0.011 

0.933081 

427S1 

1 

7 

1 15M3 

0.538 

0.933081 

427S5 

1 

6 

1 15M4 

0.130 

0.933081 

431F7 

1 

2 

115R1 

0.111 

0.933081 

431F7 

1 

2 

115R2 

0.062 

0.933081 

431R7 

1 

1 

722  . 

2833 

0  3 

121M1 

0.034 

3.913892 

325S1 

2 

20 

121M2 

0.140 

3.913892 

423S2 

2 

17 

121R1 

0.140 

3.913892 

423S2 

2 

17 

121M3 

0.044 

3.913892 

427S0 

1 

4 

121M4 

0.073 

3.913892 

427S5 

1 

6 

121M5 

0.159 

3.913892 

431F7 

1 

2 

121R2 

0.331 

3 .91  3392 

431F7 

1 

2 

121R3 

0.029 

3.913392 

431R7 

1 

1 

623.2855 


0 
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12B 

lllRl 

0.077 

4.386283 

423S0 

1 

14 

12B 

lllMl 

0.141 

4.386283 

423S0 

1 

14 

12B 

111M2 

0.008 

4 . 386283 

427S0 

1 

4 

12B 

1 1 1M3 

0.075 

4.386283 

427S5 

1 

6 

12B 

1 1 1M4 

0.178 

4.386283 

431F7 

2 

2 

12B 

1 1 1R2 

0.404 

4.386283 

431F7 

2 

2 

12B 

1 11M5 

0.041 

4.386283 

431R7 

2 

1 

12B 

1 1 1R3 

0.076 

4.386283 

431R7 

2 

1 

12C 

418  . 

9456 

0  2 

12C 

W12MM1 

0.227 

2.656809 

423S2 

2 

17 

12C 

W12MR1 

0.773 

2.656809 

423S2 

2 

17 

13A 

113. 

0546 

1  8 

13A 

131D 

0.8 

431F7 

3 

2 

13A 

131M1 

0.235 

2.219961 

423S0 

1 

14 

13A 

131R1 

0.130 

2.219961 

423SO 

1 

14 

13A 

131M2 

0.289 

2.219961 

423S4 

1 

15 

13A 

131R2 

0.141 

2.219961 

423S4 

1 

15 

13A 

13  1M3 

0.075 

2.219961 

431F7 

1 

2 

13A 

131R3 

0.069 

2.219961 

431F7 

1 

2 

13A 

131M4 

0.038 

2.219961 

431R7 

2 

1 

13A 

131R4 

0.023 

2.219961 

431F7 

2 

2 

13B 

234  . 

0806 

1  6 

13  B 

132D 

0 . 3 

431F7 

3 

2 

13B 

1  32M1 

0.177 

1  .6  14848 

423S0 

1 

14 

13B 

1  32R1 

0.110 

1  .6  14348 

423S0 

1 

14 

13B 

1  32M2 

0.269 

1  .6  14848 

423S4 

2 

15 

13B 

1  32R2 

0.103 

1  .614848 

423S4 

2 

15 

13B 

13  2M3 

0.115 

1.614848 

431R7 

2 

1 

13B 

1  32R3 

0.221 

1  .  614843 

431R7 

2 

1 

13C 

912. 

4177 

1  4 

13C 

1  34D 

0 . 8 

431  F7 

3 

2 

13C 

134M1 

0.579 

2.401157 

423S4 

2 

15 

13C 

134R1 

0.146 

2.401157 

423S4 

2 

15 

13C 

1  34H2 

0.010 

2.401157 

431  F7 

1 

2 

13C 

134R2 

0.265 

2.401157 

431  F7 

1 

2 

13D 

2164. 

4912 

1  6 

13D  133D 


0 . 8 


431F7  3  2 


133M1  0.110 
133M2  0.125 
133M3  0.082 
133R1  0.320 
133M4  0.157 
133R2  0.206 

2164 . 5000 

136D 

136M1  0.134 
136R1  0.075 
136M2  0.443 
136R2  0.254 
136M3  0.094 

1237.0099 

135D 

135M1  0.696 
135R1  0.163 
135M2  0.112 
135R2  0.029 

2134.4690 

131D 

131M1  0.235 
131R1  0.130 
131M2  0.289 
131R2  0.141 
131M3  0.075 
131R3  0.069 
131M4  0.038 
131R4  0.023 

2840.9060 

141M1  0.069 
141R1  0.121 
141M2  0.190 
141M3  0.117 
141R2  0.193 
141M4  0.212 
141R3  0.098 

853.9690 

141M1  0.069 


2 . 183050 
2.183050 
2 .183050 
2.183050 
2 .183050 
2.183050 

1  5 

0.8 

1 .528385 
1.528385 
1 .528385 
1.528385 
1 .528335 

1  4 

0.8 

2.907170 

2.907170 

2.907170 

2.907170 

1  8 

0.8 

1.540690 
1 .540690 
1 .540690 
1 .540690 
1.540690 
1 .540690 
1.540690 
1 .540690 

0  7 

2.298255 
2.298255 
2.298255 
2.298255 
2.298255 
2.298255 
2 . 298255 

0  7 

3 .930620 


423S0  1  14 
423S4  1  15 
431F7  2  2 
431F7  2  2 
431R7  2  1 
431R7  2  1 


431F7  3  2 

423S0  1  14 
423S0  1  14 
423S4  1  15 
423S4  1  15 
431F7  2  2 


431F7  3  2 

423S4  2  15 
423S4  2  15 
431R7  2  1 
431R7  2  1 


431F7  3  2 

423S0  1  14 
423S0  1  14 
423S4  1  15 
423S4  1  15 
431F7  1  2 
431F7  1  2 
431R7  2  1 
431F7  2  2 


423S0  1  14 
423S0  1  14 
427S0  1  4 
431F7  2  2 
431F7  2  2 
431R7  2  1 
431R7  2  1 


423S0  1  14 


14B 

141R1 

0.121 

3 .9306  20 

423S0 

1 

14 

14B 

141M2 

0 .190 

3.930620 

427S0 

1 

4 

14B 

141M3 

0.117 

3.930620 

431F7 

2 

2 

14B 

141R2 

0 .193 

3.930620 

431F7 

2 

2 

14B 

141M4 

0.212 

3.930620 

431R7 

2 

1 

14B 

141R3 

0.098 

3.930620 

431R7 

2 

1 

14C 

977  . 

0434 

0  7 

14C 

142M1 

0.126 

4.374974 

423S0 

1 

14 

14C 

142M2 

0.096 

4.374974 

427S0 

1 

4 

14C 

142M3 

0.128 

4.374974 

427S5 

1 

6 

14C 

142M4 

0.092 

4.374974 

431F7 

1 

2 

14C 

142R1 

0.284 

4.374974 

431F7 

1 

2 

14C 

142M5 

0.202 

4.374974 

431R7 

2 

1 

14C 

142R2 

0.072 

4.374974 

431R7 

2 

1 

14D 

806  . 

2680 

0  11 

14D 

W14GM1 

0.006 

1.475129 

423S4 

2 

15 

14D 

W14GM2 

0.056 

1 .4751  29 

426F7 

2 

11 

14D 

W14GM3 

0.007 

1.475129 

427S0 

2 

4 

14D 

W14GM4 

0.010 

1.475129 

427S5 

2 

6 

14D 

W14GM5 

0.003 

1.475129 

427S5 

1 

6 

14D 

W14GM6 

0.468 

1 . 475129 

431R7 

2 

1 

14D 

W14GM7 

0.010 

1.475129 

431F7 

1 

2 

14D 

W14GR1 

0.026 

1 . 475129 

426F7 

2 

11 

14D 

W14GR2 

0.004 

1.475129 

427S5 

1 

6 

14D 

W14GR3 

0.394 

1.475129 

431R7 

2 

1 

14D 

W14GR4 

0.016 

1.475129 

431F7 

3 

2 

14E 

793  . 

1480 

0  8 

14E 

144M1 

0.442 

2.655312 

423S4 

2 

15 

14E 

144R1 

0.112 

2.655312 

423S4 

2 

15 

14E 

144M2 

0.039 

2.655312 

427S0 

1 

4 

14E 

144  M3 

0.051 

2.655312 

427S5 

1 

6 

14E 

144M4 

0.067 

2.655312 

431F7 

2 

2 

14E 

144R2 

0.041 

2.655312 

431F7 

2 

2 

14E 

144M5 

0.158 

2.655312 

431R7 

2 

1 

14E 

144R3 

0.090 

2.655312 

431F7 

2 

2 

14F 

1117. 

3239 

0  7 

14F  145M1  0.078  2.416026  423S0  1  14 


14F 

145M2 

0.197 

2.416026 

423S4 

1 

15 

14F 

145R1 

0.096 

2.416026 

423S4 

1 

15 

14F 

145M3 

0.160 

2.416026 

427S0 

1 

4 

14F 

145M4 

0.072 

2.416026 

427S5 

1 

6 

14F 

145M5 

0.279 

2.416026 

431F7 

2 

2 

14F 

145R2 

0.118 

2 .416026 

431F7 

2 

2 

23A 


359 . 2000 


0 


7 


23  A 

W23AM1 

0.013 

2  .  7  1  2000 

325S1 

2 

20 

23A 

W23AM2 

0 .309 

2.712000 

426S7 

2 

8 

23A 

W23AM3 

0.036 

2.712000 

427S0 

1 

4 

23A 

W23AM4 

0.132 

2 .712000 

427S5 

2 

6 

23A 

W23AM5 

0.055 

2 . 7 12000 

431F7 

2 

2 

23A 

W23AR1 

0.273 

2.712000 

426S7 

3 

8 

23  A 

W23AR2 

0.127 

2.712000 

431F7 

2 

2 

23N 

1628  . 

6  6  46 

0  5 

23N 

23JM1 

0.265 

1  .  38  2284 

423S0 

1 

14 

23N 

23JR1 

0.049 

1  .  382284 

423S0 

1 

14 

23N 

23JM2 

0.278 

1  .  382284 

426F7 

2 

1 1 

23N 

23JR2 

0.238 

1  .  382284 

426F7 

2 

1 1 

23N 

23JM3 

0.170 

1.382284 

427S0 

1 

4 

23P 

3  246  . 

7517 

0  4 

23P 

23HM1 

0.484 

1  .681402 

426F7 

3 

1 1 

23P 

23HR1 

0.381 

1 .681402 

426F7 

3 

1 1 

23P 

23HM2 

0.041 

1  .68  1402 

431F7 

1 

2 

23P 

23HR2 

0.094 

1.681402 

431F7 

1 

2 

23Q 

1  329  . 

8004 

0  2 

23Q 

23GM1 

0.285 

4 .1  36834 

426F7 

3 

11 

23Q 

23GR1 

0.715 

4  .  1  36884 

426F7 

3 

1  1 

23R 

346  . 

j588 

0  4 

23R 

2  3MM1 

0.312 

3 .534767 

325S1 

2 

20 

23R 

23MR1 

0.615 

3  .  534767 

325S1 

2 

20 

23R 

23MM2 

0.019 

3.534767 

427S0 

1 

4 

23R 

23MR2 

0.054 

3  .  53476  7 

431  F7 

2 

2 

23S 

8271  . 

2920 

0  8 

23S 

W23BM1 

0.306 

2.059589 

426S7 

3 

8 

23S 

W23BM2 

0.083 

2.059589 

427S0 

1 

4 

23S 

W23BM3 

0.222 

2.059589 

427S5 

1 

6 

23S 

W23BM4 

0.028 

2.059589 

431F7 

1 

2 

23S 

W23BR1 

0.028 

2.059589 

423S0 

3 

14 

23S 

W23BR2 

0.194 

2.059589 

426S7 

3 

8 

23S 

W23BR3 

0.083 

2.059589 

431R7 

2 

1 

23S 

W23BR4 

0.056 

2.059589 

431F7 

1 

2 

23T 

305  . 

9016 

0  6 

23T 

W23LM1 

0.811 

4.594358 

426S7 

3 

8 

23T 

W23LM2 

0.010 

4  .  5943  53 

427S5 

3 

6 

23T  W23LM3  0.020  A.59A353  431 R 7  2  1 

23T  W23LM4  0.010  4.594353  43IF7  2  2 

23T  W23LR1  0.139  4.594358  426S7  3  3 

23T  U23LR2  0.010  4.594358  431R7  2  1 

23V  2238.3147  0  2 

23V  23ZM1  0.093  2.133305  426F7  2  11 
23V  23ZR1  0.907  2.138305  426F7  2  11 

41A  347.4592  0  5 

41A  W41AM1  0.394  1.884295  426F7  2  11 

41A  W41AM2  0.031  1.884295  431F7  1  2 

41A  W41AR1  0.031  1.384295  423S1  1  16 

41A  W41AR2  0.504  1.884295  426F7  2  11 

41A  W41AR3  0.040  1.884295  431F7  1  2 

41 B  1519.7465  0  11 

41B  W41BM1  0.009  3.095565  423S1  2  16 

41B  W41EM2  0.003  3.095565  423S1  3  16 

41B  W41BM3  0.245  3.095565  426F7  2  11 

41B  W41BM4  0.003  3.095565  427S5  2  6 

41B  W41BM5  0.023  3.095565  431F7  2  2 

41B  W41BR1  0.023  3.095565  423S1  2  16 

41B  W41ER2  0.634  3.095565  426F7  2  11 

41B  W41BR3  0.003  3.095565  426F7  2  11 

41B  W41BR4  0.003  3.095565  427S0  2  4 

41B  W41BR5  0.003  3.095565  427S5  1  6 

41E  W41BR6  0.051  3.095565  431F7  1  2 

41C  1736.1019  0  6 

41 C  U41CM1  0.015  3.778195  423S1  2  16 

41C  y41C.V2  0.338  3.7  73  1  9  5  426  F7  2  11 

41C  W41CM3  0.015  3.778195  431F7  2  2 

41 C  W41CR1  0.015  3.773195  423S1  2  16 

41C  W41CR2  0.537  3.773195  426F7  2  11 

41C  W41CR3  0.030  3.778195  431F7  1  2 

41D  671.2315  0  11 

41D  W41BM1  0.009  2.992214  423S1  2  16 
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Experimental  Design  And  Output 


This  appendix  contains  the  design  levels  for  the  central 
composite  design  as  well  as  the  simulation  output  obtained  by 
setting  the  factors  at  the  given  levels. 
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BMDP9R  Outpu  t 


This  appendix  Includes  BMDP9R  output  indicating  the 
'Best'  subset  of  variables  chosen  for  aircraft  availability 
(Table  D.l)  and  sortie  generation  rate  (Table  D.2). 


Table  D.l 

BMDP9R  Output  Fof  Aircraft  Availability 
STATISTICS  FOR  'BEST'  SUBSET 


MALLOWS'  CP  2.92 
SQUARED  MULTIPLE  CORRELATION  .80720 
MULTIPLE  CORRELATION  .89845 
ADJUSTED  SQUARED  MULT.  CORR .  .77754 
RESIDUAL  MEAN  SQUARE  .001226 
STANDARD  ERROR  OF  EST.  .035010 
F-STATISTIC  27.21 
NUMERATOR  DEGREES  OF  FREEDOM  2 
DENOMINATOR  DEGREES  OF  FREEDOM  13 
SIGNIFICANCE  (TAIL  PROB.  )  .0000 


NOTE  THAT  THE  ABOVE  F-STATISTIC  AND 
ASSOCIATED  SIGNIFICANCE  TEND  TO  BE 
LIBERAL  WHENEVER  A  SUBSET  OF  VARIABLES 
IS  SELECTED  BY  THE  CP  OR  ADJUSTED 


R-SQUARED 

CRITERIA. 

VARIABLE 

REGRESSION 

STAND. 

STAND. 

T- 

TOL- 

CONTRI 

BUTION 

NO.  NAME 

COEFFICIENT 

ERROR 

COEF. 

STAT. 

FRANCE 

TO  R-SQ 

INTERCEPT 
2  b 

.917869 

.0730662 

.0107198 

.0123788 

12.366 

.719 

85.62 

5.90 

1 .00000 

.51669 

14  bsq 

-  .0547960 

.0123806 

-.539 

-4.43 

1 .00000 

.29051 

THE  CONTRIBUTION  TO  R-SQUARED  FOR  EACH  VARIABLE  IS  THE  AMOUNT 
BY  WHICH  R-SQUARED  WOULD  BE  REDUCED  IF  THAT  VARIABLE  WERE 
REMOVED  FROM  THE  REGRESSION  EQUATION. 


Table  D.2 


BMDP9R  Output  For  Sortie  Generation  Rate 
STATISTICS  FOR  'BEST'  SUBSET 


MALLOWS'  CP 

SQUARED  MULTIPLE  CORRELATION 
MULTIPLE  CORRELATION 
ADJUSTED  SQUARED  MULT.  CORR. 
RESIDUAL  MEAN  SQUARE 
STANDARD  ERROR  OF  EST. 
F-STATISTIC 

NUMERATOR  DEGREES  OF  FREEDOM 
DENOMINATOR  DEGREES  OF  FREEDOM 
SIGNIFICANCE  (TAIL  PROB.) 


1  .04 
.50907 
.71349 
.43354 
.001596 
.039944 
6.74 
2 
13 
.0098 


NOTE  THAT  THE  ABOVE  F-STATISTIC  AND 
ASSOCIATED  SIGNIFICANCE  TEND  TO  BE 
LIBERAL  WHENEVER  A  SUBSET  OF  VARIABLES 
IS  SELECTED  BY  THE  CP  OR  ADJUSTED 
R-SQUARED  CRITERIA. 


VARIABLE 
NO.  NAME 


REGRESSION  STAND. 
COEFFICIENT  ERROR 


STAND. 

COEF. 


T- 

STAT. 


TOL- 


CONTRI- 

BUTION 


ERANCE  TO  R-SQ 


INTERCEPT  2.92553  .0122307  55.123  239.20 

2  b  .0383215  .0141235  .527  2.71  1.00000  .27802 

14  bsq  -.0349404  .0141257  -.481  -2.47  1.00000  .23105 


THE  CONTRIBUTION  TO  R-SQUARED  FOR  EACH  VARIABLE  IS  THE  AMOUNT 
BY  WHICH  R-SQUARED  WOULD  BE  REDUCED  IF  THAT  VARIABLE  WERE 
REMOVED  FROM  THE  REGRESSION  EQUATION. 
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